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SUMMARY 

This  report  documents  the  findings  of  a  year  long  study  to  asses  the  impact  of  Large  Woody 
Debris  (I.WD)  upon  channel  evolution  and  morphology  in  unstable  sand  bed  rivers  in  northern 
Mississippi.  The  aim  of  this  research  is  to  gain  an  improved  understanding  of  the  basin-wide 
impact  of  LWD  dynamics  in  unstable  and  stable  channel  environments  and  to  develop  a  set  of 
coherent  debris  management  strategies  for  erosion  control,  habitat  enhancement,  and 
maintenance/design  considerations  for  run-of-river  structures,  based  upon  sound  geomorphic 
and  engineering  analysis. 

Data  from  the  US  Army  Corps  of  F.ngineers  Demonstration  Erosion  Control  (DEC)  survey 
program,  conducted  in  May  1995,  has  been  used  to  locate  significant  debris  jams  with  respect 
to  planform  and  long  profile  data  on  23  river  reaches  in  northern  Mississippi  The  reaches 
surveyed  are  between  4000  and  12000  feet  long  and  range  in  upstream  basin  area  from  3.5  to 
150  square  miles.  A  comprehensive  understanding  of  debris  dynamics  can  be  attained  from 
surveying  these  channels  because  reaches  fall  into  several  categories  including,  stable/unstable 
reaches,  straight/meandering  reaches  and  reaches  which  have  either  a  predominantly 
agricultural  or  wooded  riparian  zone.  The  debris  jams  in  each  reach  have  been  surveyed  in 
detail  to  determine  the  mechanisms  and  locations  of  debris  input,  jam  impact  upon  channel 
morphology  and  sediment  routing  and  jam  stability  over  time  The  final  of  these  objectives  has 
been  assessed  by  comparing  the  survey  results  of  the  current  study  with  those  obtained  in  the 
1994/1995  research  effort  (see  Wallerstein.  1995) 

An  up  to  dale  review  of  the  literature  concerning  the  geomorphic  impacts  of  in-channel 
LWD,  and  current  LWD  management  strategies  is  presented. 

Survey  nd  reconnaissance  results  are  presented  and  analysed  Findings  show  that  the 
locations  of  debris  input  may  be  predicted  using  simple  geomorphic  variables  The  frequency 
of  jams  and  volume  of  debris  appear  to  only  be  very  weakly  related  to  drainage  basin  a'oa, 
composite  channel  width  and  unit  stream  power,  which  are  three  potentially  p  cdictive 
independent  variables.  The  distribution  of  sedimentation  and  scour  associated  with  debris  jams 
appears  to  have  an  explainable  distribution  when  related  to  drainage  basin  area  Comparison  of 
May  1994  with  May  1995  long  profile  data  sets  shows  that  the  majority  of  jams  have  remained 
in  place  over  the  intervening  period,  only  minor  jams  appear  to  have  been  displaced  and 
several  new  jams  have  been  recorded  and  largely  attributed  to  debiis  input  caused  by  bank 
instability 
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An  updated  version  of  the  LWD  Management  Program  presented  in  Wallerslein  (1995)  is 
enclosed,  on  a  disk,  with  this  document,  A  second  program  is  presented  here  (also  enclosed 
on  the  disk)  which  calculates  the  probability  of  debris  build-up  at  bridge  piers,  and  the 
associated  debris  induced  scour,  based  upon  modified  theoretical  equations  published  by 
Melville  and  Dongol  (1992)  and  Simons  and  Li  (1979), 

Preliminary  conclusions  and  management  recommendations  are  made,  based  upon  the  findings 
obtained  thus  far. 


UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  rej;t  -ri  can  be  converted  as  fol'ows 
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1  INTRODUCTION 

There  has  been  increasing  interest  in  the  role  of  vegetation  in  fluvial  geomorphology  in  recent 
years  because  it  has  been  recognised  that  river  dynamics  cannot  be  fully  understood  without 
taking  into  account  the  impact  that  vegetation  has  upon  bank  stability,  flow  velocity,  and 
riverine  habitat. 

As  a  consequence  the  study  of  in-channel  Large  Woody  Debris  (LWD)  or  Coarse  Woody 
Debris  (CWD)  as  it  is  sometimes  referred  to  (that  is  trees,  branches  and  other  larger  organic 
matter  operationally  defined  as  material  with  a  length  greater  than  1  metre)  and  its 
accumulation  as  jams  or  dams  and  impact  upon  the  channel  environment  has  become  a  topic 
receiving  increasing  research  interest  over  the  past  5  to  10  years. 

In  a  review  of  relevant  literature  undertaken  in  an  earlier  study  (Wallerstein  1994)  it  was 
established  that  a  large  proportion  of  the  research  performed  to  date  has  been  carried  out  in 
upland  areas,  and  in  stable,  gravel  bed  rivers  such  as  in  the  Pacific  North  West  (Hogan  et.  al., 
1995;  Fetherston  et.  al.,  1995)  to  determine  the  impact  of  LWD  on  salmon  habitat  and 
migration,  and  in  relation  to  logging  operations  and  forest  management.  Very  little  is  known 
about  the  impact  of  LWD  in  sand  bed,  or  unstable  rivers.  Much  of  the  work  is  fairly  qualitative 
and  observational  in  nature  and  there  has  been  little  emphasis  on  determining  the  key  variables 
at  play  in  LWD  dynamics,  and  the  modes  of  their  interaction.  Most  studies  have  also  been 
undertaken  in  isolated  reaches,  rather  than  covering  basin-wide  debris  processes,  although 
there  are  one  or  two  notable  exceptions  (see  Gregory  et  al.,  1 993). 

At  present,  LWD  management  is,  therefore,  conducted  from  an  incomplete  understanding  of 
debris  impacts  and  dynamics  and  operational  maintenance  of  debris  is  carried  out  on  an  ad  hoc 
basis. 

The  ongoing  aim  of  this  esearch  effort  is  to  asses  the  catchment  wide  impact  of  LWD  over  a 
wide  range  of  channel  sizes  but  in  unstable,  rapidly  evolving  rivers  with  sand,  clay  and  loess 
bed  and  banks.  The  research  in  this  project  lias  been  centred  on  streams  in  the  DF.C 
(Demonstration  Hiusion  Control)  watersheds  draining  the  Blutl  Line  hills  of  Northern 
Mississippi,  which  are  known  to  be  evolving  rapidly  in  response  to  complex  response  in  the 
fluvial  system  following  catchment  land-use  changes  and  past  engineering  interventions 
The  specific  aims  of  this  research  ar  c: 
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1)  To  collect  a  large,  meaningful  data  set  concerning  the  reach  scale  and  basin-wide  intlucnces 
of  LV/D  on  channel  morphology  in  a  different  type  of  channel  environment  to  that  which  has 
been  studied  so  far,  namely  unstable,  rapidly  evolving  sand-bed  rivers. 

2)  To  asses  whether  there  are  preferential  sites  of  debris  input  and  accumulation  within  the 
channel  environment  and  the  stability  of  debris  jams  in  term  of  there  longevity  in  a  particular 
reach. 

3)  To  investigate  how  effectively  debris  jams  inhibit  or  promote  bed  scour,  sediment  transport 
and  storage  in  order  to  determine  wnether  they  are  net  stabilising  or  destabilising  elements  in 
the  system 

4) To  asses  the  impact  of  debris  at  ain-of-river  hydraulic  structures  such  as  grade  controls, 
bridges,  bendway  weirs,  locks  and  dam  sluices. 

5)  To  develop  a  set  of  guidelines  for  in-channc!  LWD  management  that  could  be  used  by 
engineers  and  river  managers  as  an  aid  to  assessment,  design  and  maintenance  of  stable 
channels,  and  guidelines  for  LWI)  management  technologies  at  run-of-river  structures. 

This  report  presents  an  up  to  date  review  of  literature  concerning  the  geomorphologica! 
impact  of  LWD  and  LWD  management  strategies. 

Data  collection  was  undertaken  during  a  three  week  survey  program  in  May  IW5  with  the 
assistance  of  the  Colorado  State  University  DRC  survey  crew.  A  comprehensive 
reconnaissance  survey  was  made  of  all  debris  jams  in  the  23  study  reaches,  and  each  site  was 
also  surveyed  into  the  long  profile  and  cross  section  data  to  enable  comparison  with  the  data 
obtained  in  May  I(td4  Analysis  of  these  data  are  include  within  this  report 
The  geomorphological  characteristics  of  jams  in  each  reach  have  been  analysed  and  plotted 
against  independent  catchment  variables,  including  drainage  basin  area,  stream  power  and 
average  channel  top  width  to  determine  whether  the  geomorphological  effects  of  LWD  have  a 
coherent  and  predictable,  spatial  relationship.  Debris  jam  sediment  budgets  have  also  been 
calculated  and  related  to  spatial  parameters  to  determine  whether  the  net  impact  of  debris  jams 
is  sediment  retention  or  sediment  scour  and  mobilisation  An  understanding  of  this  factor  is 
important  as  it  will  indicate  whether  LWD  is  a  net  stabilising  or  destabilising  agent  in  sand- 
bed  livers 

This  report  also  contains  an  updated  version  of  the  LWD  Management  Program  that  was  first 
presented  in  Wallenstein  (I'tUM 


A  second  program  has  also  been  developed,  and  is  presented  here,  which  calculates  the 
probability  of  debris  build  up  at  bridge  piers,  and  the  associated  debris  induced  scour,  based 
upon  modified  theoieticai  equations  published  by  Melville  and  Dongol  (1992)  and  Simons  and 
Li  (1979)  These  two  programs  are  discussed  in  chapter  5.  The  two  programs  arc  included  on 
a  disk  inside  the  back  cover  of  this  report,  ar  i  user  manuals  are  included  in  Appendices  1  and 
2. 

The  long-term  aim  of  this  research  is  an  improved  understanding  of  the  basin-wide  impact  of 
LWD  dynamics  in  unstable  and  stable  channel  environments  and  the  development  of  coherent 
basin-wide  debris  management  strategies  for  erosion  control,  habitat  enhancement,  and 
rnainienance'design  procedure  for  DEC  and  run-of-river  structures,  based  upon  sound 
gsomomhic  and  engineering  analysis. 


2  LITERATURE  REVIEW 

2.1  INTRODUCTION 

Organic  or  woody  debris  is  an  important  channel  independent  variable  in  many  fluvial  systems 
(Hogan,  1987).  For  example,  Bevan  (1948;  quoted  in  Keller  and  Macdonald,  1995)  concluded 
that  in  the  Middle  Fork  Willamette  River,  Oregon,  woody  debris  was  responsible  for  more 
channel  changes  than  any  other  factor. 

In  a  literature  review  of  published  material  then  available,  Hickin  (1984)  suggested  that 
vegetation  may  influence  channel  processes  through  five  mechanisms 

a)  Flow  resistance 

b)  Bank  strength 

c)  Bar  sedimentation 

d)  Formation  of  log  jams 

e)  Concave-bank  bench  deposits 

Hogan  also  identified  that  the  literature  concerning  this  subject  was  of  two  main  types:  that 
dealing  with  the  indirect  influence  relations  between  vegetation,  water,  sediment  yields  and 
river  morphology,  and  that  dealing  with  the  direct  impacts  of  channel  vegetation  on  channel 
morphology 

Since  the  1980’s  the  number  of  papers  dealing  with  vegetation  in  rivers  has  increased 
markedly,  however,  including  a  number  of  studies  concerning  Coarse  Woody  Debris  (CWD), 
(Nakamura  &  Swanson,  1993),  Large  Organic  Debris  (LOD)  (Hogan,  1987)  or  Large  Woody 
Debris  (LWD),  (Smith  &  Shields,  1902)  and  its  accumulation  as  jams  or  dams  in  river 
channels 

Studies  can  be  grouped  by  topic  into  those  dealing  primarily  with 

a)  Input  processes,  distribution  and  residence  time  of  LWD 

b)  Geomorphic  significance  of  LWD 

c)  Ecological  impact  of  LWD 

The  physical  processes  involved  in  each  topic  vary  depending  upon  the  si/.e  of  the  stream 
relative  to  that  of  the  CWD  (Nakamura  vt  at,  1993) 

Most  studies  have  been  carried  out  in  essentially  stable  channel  environments  in  the  1JS  and 
Canadian  Pacific  Northwest,  the  UK,  and  New  Zealand  Instability,  in  the  form  of  landsliding, 
is  cited  by  Pearce  <$:  Watson  (1981)  as  a  means  for  debris  to  enter  channels,  but,  more 
generally,  the  study  of  debris  impacts  in  inherently  unstable  channels  has  no!  been  addressed 
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2.2  INPUT  PROCESSES,  FORMATION  AND  RESIDENCE  TIME  OF  LWD 
2.2.1  Input  Processes 

Large  Organic  Debris  enters  river  systems  by  two  main  processes;  either  from  outside  the 
channel  due  to  bank  erosion,  mass  wasting,  windthrow,  collapse  of  trees  due  ro  ice  loading  or 
biological  factors  such  as  deaih  and  litter  fall  (Keller,  1979),  or  from  inside  the  channel, 
through  erosion  and  flotation  of  emergent  and  riparian  trees  (Hogan,  1987),  (Figure  2.1). 
Fetherston  et  al.  (1995)  suggest  that  debris  inputs  are  either  “chronic  or  episodic”.  Chronic 
inputs  are  frequent  but  small  in  magnitude  and  occur  due  to  tree  mortality  and  bank  failure, 
while  episodic  inputs  are  infrequent  but  provide  a  large  amount  of  material.  Episodic  input 
processes  include  windthrow,  ice  storm,  fire  and  flood  events.  The  dominance  of  different 
input  processes  varies  widely.  For  example  45  percent  of  inpui  s  due  to  windthrow  in  the 
Lymington  Basin,  UK  (Gregory  et.  al,  1993),  while  massive  inputs  from  landsliding  of  debris 
in  a  mountain  catchment  are  reported  by  Pearce  &  Watson  (1983),  and  by  lansliding  as  a 
consequence  of  logging  operations  in  the  Queen  Charlotte  Island,  British  Columbia  by  Hogan 
et  al.  (1995).  Keller  et  al.  (1979)  suggest  that  in  low  gradient,  meandering  streams  inputs  are 
predominantly  the  result  of  bank  erosion  and  mass  bank  wasting,  windthrow  and  ice  loading, 
while  in  mountain  streams  the  main  process  is  debris  avalanche  Diehl  &  Bryan  (1994)  found 
the  dominant  input  process  to  be  bank  erosion  in  unstable  rivers  in  Tennessee  and  noted  that 
channel  instability  could  be  a  good  indicator  of  in-channel  debris  abundance.  LWD  that  has 
been  input  by  bank  erosion  can  be  identified  and  distinguished  from  that  which  has  entered  by 
other  processes  because  the  trees  will  usually  have  an  asymmetrical  root  mass  due  to 
progressive  slipping  of  the  tree  from  the  bank  into  the  channel  (Diehl  &  Bryan,  1994).  Smith 
el  al.,  ( 1993)  found  debris  input  to  be  spatially  random  However,  the  locations  of  zones  from 
which  LWD  is  supplied  will  vary  as  a  function  of  the  distribution  of  riparian  vegetation, 
streamside  topography,  channel  characteristics  and  the  prevailing  wind  strength  and  direction, 
(Fetherston  et  al.,  5 995)  It  may  therefore  be  possible  to  determine  which  are  the  dominant 
input  factors  based  on  observations  of  these  factors  and,  thereby,  predict  the  distribution  of 
major  source  areas  within  the  catchment 
2.2.2  Formation  of  Jains 

Once  in  a  channel,  debris  may  for  in  into  jams  or  dams  Jams  usually  form  ar  ound  "key  coarse 
woody  deb  is"  (Nakamura,  1903),  which  arc  usually  large,  whole  trees  that  have  entered  the 
channel  by  one  of  the  mechanisms  mentioned  above  and  which  may  he  anchored  to  the  bed  or 


banks  at  one  or  both  ends  Smaller  debris  floating  down  the  channel  then  accumulates  against 
the  key  elements,  which  acts  as  a  sieve  to  debris  J,  later  to  sediment.  If  there  is  no  fine 
debris  in  the  stream  a  mature  jam  n.ay  i  ever  form,  so  that  the  impact  of  key-debris  is  minimal. 
The  location  of  debris  jams  within  the  channel,  their  size  and  their  coherence  vary  as  functions 
of  position  in  the  catchment.  In  small  streams  much  debris  will  accumulate  where  it  falls 
because  the  flow  is  not  competent  to  move  coarse  material,  and  it  is  in  larger  streams  that 
distinct  jams  may  form.  In  yet  larger  rivers  debris  may  never  accumulate  because  it  is  carried 
away  downstream.  Piegay  (1993)  observed  debris  distribution  in  a  sixth  order  river  in  France 
and  found  that  most  material  was  deposited  on  the  channel  margins,  forming  a  narrow  debris 
line  rather  than  in-flow  jams  Wallance  &  Benke  (1984)  noted  a  similar  distribution  in 
meandering  rivers  in  the  south  east  USA  where  dense,  partial  jams  formed  at  a  angle  to  the 
main  flow.  As  mean  channel  dimensions  and  flow  competence  increase  downstream  more  and 
more  debris  will  be  moved  from  its  position  of  input,  until  all  but  the  hugest  trees  are 
fransported  This  process  relationship  may  r  suit  in  a  trend  of  reducing  LWD  frequency 
downstream,  but,  at  the  same  time,  an  increase  in  the  volumetric  size  of  each  jam  (Swanson  et 
al.,  1982). 

2.2.3  Residence  time  of  debris  jams 

The  residence  time,  or  persistence,  of  debris  jams  is  an  important  factor,  which  determines  the 
timespan  over  which  channel  morphology  at  a  jam  site  will  be  affected.  The  influence  exerted 
by  jams  on  channel  morphology  also  varies  with  time  as  the  debris  in  the  jam  structure 
deteriorates  (Hogan  et.  al.,  1995).  Assessing  residence  time  is  difficult  and  estimates  range 
between  12  months,  for  a  36%  change  or  removal  (Gregory  &  Gurnell  (1985),  to  40-90  years 
(Hogan,  1987),  to  200  years  for  streams  in  British  Columbia  (Keller  &  Tally,  1979)  Residence 
times  may  vary  as  a  function  of  drainage  basin  area,  and  are  largely  dependent  upon  the  return 
period  of  a  flood  with  a  magnitude  which  is  capable  of  entraining  a  significant  proportion  of 
the  trapped  debris  or  moving  larger  key  components  of  the  jam.  Other  important  factors 
affecting  jam  persistence  are  average  tree  dimensions  and  wood  deterioration  rate.  Swanson  et 
al.  (1982)  discovered  that  the  density  and  volume  of  in-channel  debris  are  greater  in  rivers 
which  flow  through  coniferous  forests,  than  it  is  in  those  that  flow  though  deciduous  forests 
This  is  because  conifers  arc,  on  average,  taller  and  have  slower  decay  rates  than  deciduous 
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Figure  2.2  Energy  transformations  in  streams  with  log  steps  (adapted  from  Marston, 
1982)  ,  , 


Potential  stream  energy  per  unit  mass  of  water  (PE/m)  is  directly 
proportional  to  h,  or  the  relief  in  a  specific  stream. 

PE  dissipation  by  log  steps  -  Cumulative  change  of  water  suiface 
elevation  (hi  +h2+  etc)  as  a  percentage  of  total  stream  relief  (0 
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Log  steps  cause  a  reduction  in  potential  energy  that  would 
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„  2.3  GEOMQRPHIC  SIGNIFICANCE  OF  LWD 

2.3.1  fleets  of  channel  scale 

It  is/mportant  to  recognise  that  processes  are  scale  dependent  and  that  the  influence  of  LWD 
or''  channel  and  valley  morphology  may  change  systematically  downstream  through  the 
network  (Abbe  &  Montgomery,  1993).  Zimmerm.n  et  al  (1967)  found  that  debris 
accumulations  in  a  very  small  stream  completely  obscured  the  usual  hydraulic  geometry 
relations,  while  Robinson  &  Boschta  (1990),  and  Keller  &  Tally  (1‘  79)  suggest  that  debris 
loadings  increase  with  stream  size.  Gregory  et  al  (1985),  have  characterised  jams  into  three 
types 

1 )  Active  (form  a  complete  barrier  to  water  and  sediment  movement,  and  create  a 

distinct  step  or  fall  in  the  channel  profile) 

2)  Complete  (a  complete  barrier  to  water/sediment  movement,  but  no  step  formed) 

3)  Partial  (only  a  partial  barrier  to  flow) 

They  suggest  that  these  types  become  sequentially  more  prevalent  as  channel  size  ini  l  eases. 
Once  trees  fall  into  a  stream,  their  influence  on  channel  form  and  process  may  be  quite 
different  from  that  when  they  were  on  the  banks,  changing  from  tabilising  to  destabilising 
through  causing  local  bed  scour  and  basal  erosion  of  the  banks  Thus,  jams  represent  a  type  of 
auto-diversion,  that  is,  a  change  in  channel  morphology  triggered  by  the  fluvial  process  itself 
(Keller  &.  Swanson,  1979).  The  type  and  degree  of  impact  on  channel  morphology  depends 
primarily  on  the  channel  width/trec  height  ratio  and  on  debris  orientation  relative  to  the  flow. 
Mean  discharge  and  the  dominant  discharge  recurrence  interval  are  also  important  because  the 
higher  the  flow  is  relative  to  jam  size,  the  smaller  will  be  the  jam’s  impact  in  terms  of  acting  as 
a  How  diverter  and  roughness  element  The  principal  effects  of  debris  upon  channel 
morphology  arc  described  below 

2.3.2  Impact  of  debris  jams  upon  channel  morphology 

LWD  influences  the  geomorphology  of  rivers  on  three  levels  (Gray,  1974),  the  overall  channel 
form;  detailed  features  of  the  channel  topography;  and  channel  roughness 
Heede  (1985),  Smith  (1993),  Andrus  et  al.  (1988)  and  Mosley  (1981)  have  all  observed  that 
the  spatial  distribution  and  number  of  pools,  riffles  and  gravel  bars  is  positively  related  to  the 
distribution  and  volume  of  LWD  in  the  channel  This  relationship  has  been  explained  through 
laboratory  experiments  by  Smith  &  Beschta  ( 1 994),  who  found  that  the  pool-riflle  sequence  in 


Figure  2.3  Schematic  diagrams  of  pool  types  (Modified  from  Robinson  and  Beschta, 
1990) 
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gravel-bed  rivers  is  maintained  by  a  combination  of  mean  boundary  shear  stress  and 
intermittent  lit)  and  drag  forces  due  to  velocity  fluctuations  aiound  debris  Random  debris 
input  will  also  distort  the  pool-riffle  sequence,  making  it  less  systematic,  so  that  the  long- 
profilc  has  very  little  spatial  memory,  or  periodicity  (Robinson  &  Beschta,  1990)  Robinson 
and  Beschta  ( 1990)  devised  a  pool  classification  system,  containing  six  pool  categories  (lateral 
scour,  fluvial,  plunge,  underflow,  d.  Hector  and  dam)  based  on  How  and  debris  interaction  (see 
figure  2..1).  Other  studies  have  shown  that  a  considerable  proportion  of  the  vertical  fall  of 
channels  an  occur  at  the  sites  of  debris  jams,  accounting  for  a  4%  of  the  vertical  drop  along  a 
412m  n  ;h  of  channel  in  Vermont  (Thompson,  1995)  and  <>0%  of  the  total  drop  in  Little  l.ost 


in 


Man  Creek  in  Northern  California  (Keller  &  Tally,  1079)  Debris  jams,  therefore,  act  as  local 
base  levels  and  sediment  storage  zones  which  provide  a  buffer  in  the  sediment  routing  system 
(Heede,  1985,  Bilby,  1981).  Thompson  (1995)  found  that  TWD  causes  an  important  negative 
feedback  mechanism,  where,  in  the  case  of  channel  degradation,  there  is  an  increase  in  debris 
input  due  to  mass  bank  failure,  which  in  turn  causes  greater  sediment  storage  Channel  bed 
elevation  is  consequently  raised  once  more  and  the  rate  of  bank  failure  and  debris  input  is 
thereby  reduced.  On  this  basis,  Klein  et  al.  (1987)  argue  that  jam  removal  can  reduce  the  base 
level  for  the  channel  upstream  and  may  trigger  bank  erosion.  However,  in  an  experimental 
study  by  Smith  et  al.  (1993a  and  b)  *t  was  found  that,  while  the  removal  of  debris  from  a  small 
gravel  bed  stream  initially  caused  a  four  fold  increase  in  bed  load  transport  at  bankfull  flow, 
the  associated  loss  of  scour  turbulence  and  greater  flow  resistance  imparted  by  alternate  bars 
actually  resulted  in  a  reduction  in  stream  power  which  was  compensated  for  by  sediment 
deposition  and  net  channel  aggredation. 

Potential  energy  is  dissipated  at  jams,  with  energy  loss  being  as  much  as  6%  o<’  total  potential 
energy  (MacDonald  et  al.,  1982).  Shields  &  Smith  (1992)  found  that  the  Darcy-Weisbach 
friction  factor  was  40'J  %  highei  at  base  flow  in  an  uncleared  river  reach  compared  to  a  clear 
condition,  but  that  this  value  declined  to  35%  at  high  flows  The  velocity  distribution  is  also 
far  more  heterogeneous  in  debris-filled  reaches,  especially  at  low  flow  Changes  of  stream 
power  distribution  due  to  flow  resistance  effects  in  turn  give  jams  the  ability  to  influence  the 
location  of  erosional  and  depositional  processes  Also  the  backwater  etfect  created  by  jam 
back-pools  may  induce  local  silting  (Keller  et  al.  1976).  Thus,  in  small,  stable  channels,  log 
steps  generally  increase  bank  stability  and  reduce  sediment  transport  rales  by  creating  falls, 
runs  and  hydraulic  jumps. 

Figure  2.2  shows  how  potential  energy  is  lost  through  a  log-step  sequence  as  outlined  by 
Marston  (1982).  The  localised  dissipation  of  energy  can,  I  vever,  result  in  associated  local 
scour  and  bank  erosion  which  causes  channel  widening.  Bank  failure  may  also  occur  through 
flow  diversion  around  a  debris  obstruction  (Murgatroyd  &  Ternan,  1983).  Davis  &  Gregory 
(1994)  have  also  suggested  a  mechanism  whereby  bank  failure  is  induced  through  the  erosion 
••.f  a  porous,  gravel,  bank  subsurface  due  to  the  greater  hydrostatic  pressure  caused  by  debris 
dammed  flow  Conversely,  Keller  &  Tally  (1979)  have  observed,  that  flow  convergence  under 
logs  may  cause  channel  narrowing,  with  sediment  storage  upstream  and  a  scour-pool 
downstn  am  of  the  log  step. 


As  drainage  area  increases,  and  the  channel  width/tree  size  ratio  exceeds  unity,  How  is  diverted 
laterally,  inducing  bank  erosion  through  local  basal  scour.  Hogan  (1987)  found  that  in 
undisturbed  channels  in  British  C'olumbia  organic  debris  orientated  diagonally  actoss  the 
channel  resulted  in  high  width  and  depth  variability.  However,  in  catchments  where  there  had 
been  logging  operations  the  majority  of  in-channel  discarded  timber  was  orientated  parallel  to 
the  flow  ind  it  subsequently  became  incorporated  into  the  stream  banks,  protecting  them  from 
erosion  Nakamura  &  Swanson  (199.1)  and  Keller  &  Swanson  (1979)  have  suggested  t  i 
there  is  a  progression  of  types  of  interaction  between  debris  jam  and  channel  processes, 
ranging  from  local  base  level  control  and  possible  local  widening  in  low-order  streams,  to 
lateral  channel  shifts  and  even  meander  cut-off  in  middle-order  channels,  where  debris  is 
moved  into  larger  more  coherent  jams  which  may  either  increase  or  decrease  the  channel 
stability  depending  upon  the  erodibility  of  bed  and  banks.  In  larger  channels  still,  bars  may 
form  and  flow  bifurcate  around  debris  obstructions.  This  last  process  has  been  documented  by 
Nanson  (1981)  in  British  Columbia,  who  found  that  organic  debris  deposited  at  low  flow 
provided  the  nuclei  for  development  of  scroll  bars,  through  the  local  reduction  of  stream 
power.  Hickin  (1984)  also  observed  crib-like  bar-head  features,  but  was  undecided  regarding 
whether  the  debris  caused  bar  formation,  or  whether  the  bars  pre-dated  and  trapped  the  debris. 
In  either  cast,  organic  debris  would,  at  the  very  least,  enhance  sediment  deposition  and  bar 
formation. 


2.4  ECOLOGICAL  IMPACT  OF  LWD 

LWD  dams  are  very  important  in  small  stream  ecosystems  because  they  provide  a  source  of 
organic  matter  and  retain  floating  leaves  and  twigs  in  the  dam  structure  and  backwater  pools. 
This  coarse  particulate  organic  matter  (CPOM)  is  broken  down  in  the  low  energy  pool 
environment  by  shredder  invertebrates,  creating  fine  particulate  organic  matter  (IPOM)  and 
dissolved  organic  matter  (POM),  which  are  the  required  energy  sources  of  a  succession  of 
invertebrate  species  who  are,  in  turn,  the  energy  source  of  high  fauna  species.  Bilby  &  Likens 
Vi  vow;  iwunu  nidi  mu  puiuuma^u  ui  mu  Maiming  siock  oi  organic  mailer  reiameci  oy  jams 
changed  from  75%  in  first  order,  to  58%  in  second  order  to  20%  in  third  order  streams 
because  the  prevalence  of  dam  type  jams  declined  downstream  The  volume  of  (’POM 
therefore  declines  downstream,  while  the  volume  of  IPOM  and  DOM  increases  This  gives 
rise  to  a  spatially  varied  invertebrate  community,  changing  from  shredders  in  small  channels  to 
gathers  of  I  POM  downstream  Smock  cl  al.  (19X9)  and  Wallanee  &  Benke  (1984)  found 
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similar  correlations  between  debris  volume  and  invertebrate  abundance  in  sand-bed  streams, 
where  debris  provides  the  only  stable  substrate  for  organic  matter  retention  and  invertebrate 
habitat  Higher  species,  such  as  fish,  use  debris  and  associated  pools  for  shade,  protection 
from  predators,  feeding  and  spawning  grounds.  The  pools  and  falls  created  by  log  step ,  also 
help  to  oxygenate  the  flow,  and  provide  a  variety  of  different  energy  environments  which  are 
can  be  colonised  by  niche  species 

In  addition  to  providing  essential  fauna  habitat,  LWD  is  also  a  vital  factor  in  the  development 
of  the  riparian  forest  mosaic  (Fetherston  et.  al.,  1995).  Debris  deposition  in  the  channel  and  on 
the  floodplain  creates  sites  of  low  boundary  shear-stress  where  vegetation  colonisation  can 
take  place  This  leads  to  the  development  of  vegetation  stabilised  islands  and  bars  (affecting 
the  geomorphological  development  of  the  channel)  which  may  subsequently  coalesce  and/or 
become  attached  to  the  bankline  to  form  new  areas  of  forested  floodplain  tl  *  provide  shade, 
hank  stability  and  supply  and  storage  of  organic  matter,  sediment,  water  and  new  LWD. 
Figure  2.4  shows  a  modified  version  of  the  LWD-drivcn  model  of  Fetherston  et  al  .,  (1995)  for 
riparian  forest  development,  based  upon  research  findings  from  the  Pacific  Northwest 


Figure  2.4  Conceptual  model  of  montane  riparian  forest  development.  (Modified  from 
Fetherston  et  al.,  1995) 


2.5  MANAGEMENT  STRATEGIES 

Until  basic  research  concerning  in-channel  LWD  began  to  suggest  otherwise,  it  was  commonly 
believed  that  LWD  was  detrimental  to  the  fluvial  system,  hydraulically,  ecologically  and 
geomorphically.  On  this  bas's,  reasons  for  debris  removal  included  : 

a)  To  improve  navigation, 

b)  To  increase  channel  conveyance  by  reducing  roughness; 

c)  To  eliminate  bank  erosion; 

d)  To  facilitate  the  migration  offish,  especially  salmon  (MacDonald,  1982) 

It  is  now  recognised  that  there  are  advantages  to  be  gained  by  maintaining  or  even  increasing 
in-channel  debris  accumulations  (Gregory  &  Davis,  1992;  Keller  &  McDonald,  1995). 
Management  strategies  that  are  currently  advocated  vary  widely,  however  This  perhaps 
reflects  our,  as  yet,  incomplete  understanding  of  LWD  dynamics  in  different  channel 
environments,  and  because  goals  vary  between  different  management  strategies.  In  this  respect 
effective  debris  management  depends  on  the  underlying  aims  of  the  proposed  management 
action 

Successful  management  also  depends  upon  a  comprehensive  understanding  of  the  following 
hydrogeomorphological  factors  (Gregory  &  Davis,  1992) 

a)  The  relationship  between  river  channel  processes  and  river  channel  morphology, 

b)  Awareness  of  the  timescales  over  which  river  channels  may  adjust; 

c)  Consideration  of  channel  management  in  the  wider  context  of  river  basin 


management 

More  specifically,  debris  management  must  consider  : 

a)  C  hannel  stream  power  characteristics, 

b)  Sediment  movement  and  storage  relationships  (high/low;  fine/coarse  sediment; 

suspended/bedload); 

c)  Channel  stability; 


d)  Size  and  character  of  river  ch-mnc 


debris  si /.( 


e)  Spacing  and  frequency  of  jams, 

f)  Size  and  character  of  jams,  and  orientations  of  component  material, 

7)  Age  and  stability  of  component  materials. 

In  an  evaluation  of  soft  engineering  for  instream  structures,  including  some  using  woody 
debris,  to  mitigate  the  effects  of  highway  construction  in  British  Columbia,  Miles  ( 1995)  found 


that  nearly  50%  of  the  structures  had  been  severely  damaged  after  8  to  14  years  Miles 
attributed  this  problem  to  insufficient  understanding  and  consideration  of  the  stability  of  the 
structures  in  a  high  energy  river  environment.  He  concludes  that  soft  restoration  techniques 
may  not  be  appropriate  in  highly  energetic  mountain  rivers,  and  that  if  restoration  is  to  be 
performed,  funding  must  be  made  available  for  long  term  monitoring  and  maintenance 
There  appears,  in  general,  to  be  a  consensus  of  opinion  amongst  researchers  interested  in 
LWD  regarding  appropriate  management  approaches  for  channel  restoration.  Bren  (1993)  and 
Nunnally  (1978)  argued  that  the  riparian  zone  should  be  left  undisturbed,  in  a  natural  state 
(although  defining  natural  is  difficult  in  most  channels),  and  that,  because  debris  is  so 
important  for  the  river  ecosystem,  debris  jams  should  be  left  in  place.  Keller  and  McDonald 
(1995)  studied  catchments  which  had  been  disturbed  by  logging  operations  They 
recommended  that  a  riparian  buffer  strip  should  be  left  to  maintain  the  natural  LWD  supply 
and  warned  that  landsliding  events  caused  by  badly  controlled  logging  operations,  can  cause 
excessive  LWD  input  which  is  detrimental  to  stream  habitat  and  flow  and  sediment 
conveyance.  There  may  be  a  case  in  streams  lacking  a  wooded  riparian  strip  for  the 
introduction  of  debris  jams  (Keller  &  McDonald,  1995).  If  a  debris  recharge  policy  is  to  be 
implemented,  however,  it  is  important  that  debris  jam  volume  and  orientation  emulates  the 
values  which  would  be  found  under  natural  conditions  (Robinson  &  Beschta,  1990).  Wallace 
&  Bcnke  (1984)  concluded  that,  in  most  instances,  the  best  management  is  probably  no 
management  except  where  adjacent  floodplains  have  to  be  protected  from  flooding 
Comprehensive  studies  of  coarse  woody  debris  in  relation  to  river  channel  management  have 
been  carried  out  by  Gregory  and  Davis  (1992)  and  Gurncll  and  Gregory  (1995a  and  b)  The 
collation  of  analyses  from  twenty  two  research  papers  with  primary  field  studies  carried  out  in 
the  New  f  orest,  UK,  Gregory  &  Davis  (1992)  demonstrated  the  significance  of  LWD  to 
channel  morphology,  processes  and  ecology  (figure  2.5)  and  produced  a  preliminary  table  of 
debris  management  criteria  based  upon  their  findings  (figure  2,6).  They  conclude  that  "...  a 
conservative  approach  to  debris  removal  should  be  adopted  for  most  areas,  but  that  different 
strategies  are  needed  according  to  the  characteristics  of  particular  localities"  (Gregory  and 
Davis,  1992,  pg.  113). 

It  should  be  noted,  however,  that  this  study,  in  common  with  most  others  cited,  was  carried 
out  in  an  essentially  stable,  equilibrium  channel  environment  where  changes  to  channel 
morphology  ate  negligible  and  significant  impacts  i elate  mostly  to  ecological  habitat  diversity 

l  s 


I 


> 


Also,  little  attention  is  paid  to  "different  strategies"  that  may  be  required  in  contrasting  channel 
environments  and  there  is  no  discussion  of  conflicts  between  practices  advocated  by  various 
organisations  in  the  USA.  For  example,  Gregory  &  Davis  (1992)  suggest  that,  based  on  their 
literature  survey,  no  debris  should  be  removed  from  channels  exhibiting  low  stability  (Figure 
2.7).  However,  this  contradicts  the  practice  described  by  Brookes  (1985,  pg.  64),  "In  North 
America  the  concept  of  channel  restoration  was  developed  in  North  Carolina  under  the 
funding  of  the  Water  Resources  Research  Institute  of  the  State  University  Restoration  is 
achieved  by  removing  debris  jams  and  providing  uniform  channel  cross-sections  and  gradients 
whilst  preserving  meanders,  leaving  as  many  trees  as  possible  along  the  stream  banks,  and 
stabilising  banks  with  vegetation  and  rip-rap  where  necessary 

Similar  approaches,  have  been  documented  and  carried  out  by  numerous  researchers  and 
organisations  in  the  USA,  including;  McConnel  et  al.  (1980),  based  upon  work  on  the  Wolf 
River,  Tennessee,  the  American  Fisheries  Society  (1983),  in  a  publication  entitled  “Stream 
Obstruction  Removal  Guidelines”,  (see  Figure  2.8),  Shields  and  Nunnally  (1984);  and  Palmiter 
(Institute  of  Environmental  Sciences,  1982). 

The  recommendations  of  Palmiter  (1982)  include  the  following 

a)  Removal  of  log-jam  material  by  cutting  it  to  a  manageable  size; 

b)  Protection  of  eroding  banks  using  brush  piles  and  log-jam  material,  with  rope  and 

wire, 

c)  Removal  of  sand  and  grave!  using  brush-pile  deflectors; 

d)  Revegetation  to  stabilise  banks  and  shade-out  aquatic  plants; 

e)  Removal  of  potential  obstructions  such  as  trees  and  branches; 

Willcke  &  Baldwin  ( 1984)  assessed  the  Palmiter  techniques  and  found  them  suitable  for  areas 
experiencing  chronic,  low  intensity  flooding  and  bank  erosion,  but  not  advisable  for  rivers  with 
extreme  flood  problems.  They  are  also  found  to  be  largely  ineffective  for  erosion  control 
where  the  mechanism  of  bank  failure  is  that  of  mass  wasting  rather  than  tractive  force  erosion 
(Hassol  wander,  1989) 

It  is  evident  from  the  preceding  discussion  of  UWD  management  strategies  that 
recommendations  vary  considerably,  from  limited  or  no  interference,  to  total  clearance  of 
debris  from  the  channel  These  apparently  contradictory  recommendations  must  be  viewed  in 
the  light  of  the  overall  management  programme  that  they  were  designed  for,  as  i equipments 
for  habitat  enhancement  differ  from  those  for  flood  de'encc. 


In 


Finally,  and  of  great  importance,  is  the  fact  that  the  recommendations  of  type  made  by 
Palmite  and  others,  address  debris  management  predominantly  in,  low  gradient,  sand-bed,  and 
perhaps  unstable,  flood  prone  rivers  (South  East  USA),  while  those  prescribed  by  Gregory  and 
Davis  (1992)  and  others  are  based  upon  findings  from  upland,  even  montane,  gravel-bed 
rivers  and  streams  (Pacific  Northwest  USA)  Process  relationships  between  the  debris  and  the 
channel  are  likely  to  differ  between  these  two  types  of  fluvial  environment,  although,  as  yet, 
these  differences  have  not  been  recognised  or  investigated.  Indeed,  while  there  is  a  w'ealth  of 
research  concerning  the  geomorphological  impacts  of  LWD  in  upland  gravel-bed  rivers,  there 
has  been  little  comparable  research  in  lowland,  sand-bed,  and/or  unstable  river  environments. 
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Hgure  2.6  Determinants  for  a  management  strategy  for  rivers  in  woodland  areas 
(modified  from  Gregory  and  Davis,  1992) 
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Figure  2.7  Definition  of  Stream  Obstruction  Conditions  (Modified  from  American 
Fisheries  Society,  1983) 


Condition  One  Management  Criteria 

These  stream  segments  have  acceptable  How  and  no  work  No  work  to  be  conducted. 

would  be  required.  I  lie',  may  contain  various  amounts  of 

instreum  debris  and  fine  sediment,  such  as  silt,  sand, 

gravel,  nibble,  boulders,  logs  and  hmsli.  In  certain 

situations  How  may  be  impeded,  but  due  to  stream  and  land 

classification  or  adjacent  land  use,  this  is  not  a  problem _ 


Condition  Two 

These  stream  segments  currently  have  no  major  llow 
impediments,  but  existing  conditions  arc  such  thut 
obstniclions  arc  likclv  to  form  in  the  near  future,  causing 
unacceptable  piobleins.  This  condition  is  generally 
characterised  by  small  accumulations  of  logs  und/or  other 
debris  which  occasionally  sjian  the  entire  stream  width. 
Accumulations  arc  isolated,  not  massive  and  do  not 
presently  cause  upstream  ponding  damage. 


Condition  Three 

l"hcsc  stream  segments  have  unacceptable  flow  problems. 
Obstructions  are  generally  elia  v'ensed  by  huge 
accumulations  of  lodged  trees,  root  .ads,  and/or  oilier 
debris  tliat  frequently  sjiuii  the  entire  stream  width 
Although  impeded,  some  llow  moves  through  the 
obstruction,  large  amounts  of  sediment  have  not  covered  or 
lodged  in  the  obstruction 


Condition  Four 

These  stream  segments  aie  characterised  by  majoi 
blockages  causing  unacceptable  How  problems. 
Obstructions  consist  o I  compacted  debris  amt/or  debris  thn 
severely  restricts  llow 


Condition  Five 

ihese  stream  segments  posses  unique,  sensitive,  or 
especially  valuable  biotic  resources  and  should  lx;  dealt 
with  on  a  case-by-case  basis.  Kxainplcs  include,  but  are  not 
limited  to:  Areas  huilxmimg  laic  or  cndangcicd  species, 
shellfish  beds,  lisli  spawning  and  icanng  areas,  and 
nxikerics. 


Management  Criteria 

Hquipment  that  will  cause  the  least  damage  to  tile 
environment  shall  lx;  selected  for  performing  Lire  work, 
hirst  consideration  will  be  given  to  the  use  of  hand  operated 
equipment  such  us  axes,  chain  saws,  and  winches  to  remove 
accumulations.  Moats  with  motors  may  bo  used  where 
needed.  When  the  use  of  hand  operated  tixils  is  not  feasible, 
heavier  equipment  may  be  used,  e.g.  small  tractors, 
backhocs,  bulldozers,  log  skidders  and  low  PSI  equipment. 
Hquipment  shall  x:  operated  in  a  manner  that  results  in 
least  damage  to  vegetation  und  soils  of  the  project  area  In 
some  cases  explosives  may  lx:  used  resulting  in  less 
damage  Debris  designated  for  removal  from  the  stream  or 
floodway  should  he  removed  or  secured  in  such  a  manner  as 
to  restrict  its  re-entry  into  the  channel,  (h  ncially,  it  should 
be  positioned  so  ns  to  induce  flixxl  llow  impediment _ 

Management  Criteria 

Hquipment  limitations  will  lie  the  same  us  for  condition  two 
segments.  Work  shall  lx:  accomplished  within  the  clianncl 
or  from  one  side  of  the  ..lunnel  where  possible.  Selective 
tree  clearing  shall  he  limited  (o  the  minimum  necessary  for 
equipment  access  and  efficient  operation  of  equipment  on 
tlie  wori.ed  side  ol  die  channel  Disposal  of  equipment  may 
>>■  amiplishcd  by  removing  it  from  the  lloodpluin  or  by 
bujuu  „  burying  01  piling,  as  appropriate,  with  the 
M.0'  •  ,;n  linen  n'  of  disturbance  to  vegetation  Piled  debris 
p  Ircquent  intervals  and  at  all  tributaries 
.  1st  tributaries. _ 

r.  Management  Criteria 

.  Ikickugc  removal  may  employ  any  equipment  necessary  to 
accomplish  the  work  in  the  h  ast  damaging  manner.  Work 
should  be  accomplished  from  one  side  ol  (lie  channel  where 
piuctiuil.  Material  shall  1)e  disused  in  accordance  willi 
guidelines  presented  above  lor  condition  three  segments. 
Spoil  piles  should  l>e  constructed  as  high  as  sediment 
properties  allow  The  placement  ol  spoil  around  the  base  ot 
mature  trees  should  lx-  avoided. 

Management  Criteria 

Special  provi  on  for  protecting  unique,  sensitive,  or 
productive  biotic  icsourecs  shall  lx:  developed  by 
appropriate  piofessionalx  01  a  case  bv  case  basis. 
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3  RESEARCH  METHODS 

Field  data  collection  for  the  current  project  was  carried  out  in  May/June  1995  with  the  co¬ 
operation  of  the  DEC  survey  crew  front  Colorado  State  University.  Twenty  three  channel 
reaches  were  surveyed  in  the  DEC  monitoring  catchments,  which  are  located  in  the  Yazoo 
Basin,  northern  Mississippi.  Figure  3.2  shows  a  geological  map  of  Mississippi  with  the  DEC 
project  area  marked.  Figure  3.3  shows  the  rivers  and  major  catchments  of  the  project  area  in 
more  detail. 

Reaches,  of  between  4000  and  12000  ft,  were  surveyed  on  the  following  creeks  : 

Nolehoe  Creek  Sarter  Creek  Lick  Creek 

Burney  Branch  James  Wolf  Creek  Long  Creek 

Sykes  Creek  Hotopha  Creek  Fannegusha  Creek 

Worsham  Creek  (East)  Worsham  Cre''k  (Middle)  Worsham  Creek  (West) 

Abiaca  Creek  Harland  Creek  Red  Banks  Creek 

Otoucalofa  Crec \  Coila  Creek  Lee  Creek 

Pen  /  Creek  Hickahala  Creek  Marcunt  Creek 

A  full  description  of  the  DEC  monitoring  site  characteristics  is  given  in  Watson  et  al  (1993). 
These  surveys  provided  a  comprehensive  data-set,  which  not  only  covers  a  range  of  drainage 
basin  areas,  from  4  to  100  miles  square,  but  also  allows  comparison  of  debris  loadings 
between  reaches  with  wooded  and  agricultural  riparian  zones,  between  straight  and 
meandering  reaches  and  between  highly  unstable  and  stabilising  or  equilibrium  reaches, 

Debris  jam  sites  have  been  surveyed  into  the  thalweg  and  cross  section  data  for  each  creek  so 
that  there  position  and  associated  changes  in  local  channel  geomorphology  can  be  monitored 
over  time.  Data  from  the  current  survey  has  been  processed  and  overlaid  and  compared  with 
that  collected  in  May  1994  so  that  an  assessment  can  be  made  of  the  rate  of  debris  input,  of 
the  longevity  of  jams  and  therefore  their  effectiveness  as  geomorphological  channel  controls, 
and  the  changing  patterns  of  associated  sedimentation  and  erosion 

Geomorphological  reconnaissance  was  also  carried  out  at  each  jam  site  to  document  the 
volume  of  debris  in  each  jam,  its  mode  of  input  into  the  channel,  to  determine  the  jam  type  in 
terms  of  impact  upon  How  pattern  and  erosion,  and  to  measure  the  volume  of  sediment 
retained  in  backwater  areas  or  bars.  The  following  variables  were  assessed  at  each  jam  site; 

1)  Debris  jam  volume  :  Estimated  volume  of  woody  material  (nr')  in  each  jam.  values  are  then 
summed  for  each  survey  reach 


2)  Morphological  classification  :  A  debris  classification  system,  modified  from  one  developed 
by  Robinson  and  Beschta  (1989),  which  describes  the  geomorphological  impact  of  debris  jams 
throughout  the  drainage  network.  Figure  2.3  shows  the  original  classification  scheme  and 
figure  5. 1  shows  the  modified  version  Jam  classification  types  arc  as  follows: 

Underflow  jams  :  In  small  catchments  where  fallen  trees  span  the  channel  at  bank-full  level 
Local  bed  scour  may  occur  under  debris  at  high  flows,  otherwise  the  in-channel 
geomorphic  impact  of  the  LWD  is  minimal. 

Dam  jams  :  In  channels  where  the  average  tree  height  to  channel  width  ratio  is  rough  equal  to 
one,  so  that  debris  completely  spans  the  channel  cross-section.  This  type  of  jam  causes 
significant  local  bank  erosion  and  bed  scour  due  to  flow  constriction,  and  backwater 
effects  will  cause  sediment  deposition  in  the  lower  energy  environment  upstream  Bars 
may  also  form  immediately  downstream  of  the  jam. 

Deflector  jams  :  Found  where  input  debris  does  not  quite  span  the  channel  so  that  flow  is 
deflected  against  one  or  both  of  the  banks  causing  localised  bed  scour  and  bank 
erosion.  Subsequent  bank  failure  results  in  the  input  of  new  LWD  material  to  the  reach 
so  that  the  jam  builds  up  further  Backwater  sediment  wedges  and  downstream  bars 
may  form  at  this  type  of  jam  provided  that  stream  power  is  dissipated  by  the  jam  below 
the  critical  level  for  the  bed  load  and  suspended  sediment  transport 
l'low  Parallel  jams  :  Found  where  channel  width  is  significantly  greater  than  the  key-debris 
length,  and  flows  are  competent  enough  to  rotate  oebris  so  that  it  lies  parallel  to  the 
flow  Debris  is  also  transported  downstream  in  high  (lows  and  deposited  against  the 
bank-base  on  the  outside  of  meander  bends  or  at  channel  obstructions  such  as  man¬ 
made  structures  Related  bank  erosion  and  bed  scour  will  be  minimal,  and  bank  toes 
may  even  be  stabilised  by  debris  build-up.  Flow  parallel  debris  may  also  initiate  or 
accelerate  the  formation  of  mid-channel  and  lateral  bars. 

3)  Blockage  classification:  Jams  are  classified,  using  the  scheme  developed  by  Gregory  et  al 
(1985),  according  to  their  potential  to  block  the  downstream  movement  of  water  and 
sediment.  The  classification  types  are: 

Active'.  Jam  forms  a  complete  barrier  to  water  and  sediment  movement  and  also  creates  a 
distinct  step  or  fall  in  the  channel  profile. 

<  'ample  te  complete  barrier  to  water/sediment  movement,  but  no  significant  step. 

Partial.  Jam  is  ally  a  partial  barrier  to  (low. 
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4)  Alpha/  Beta  Indices:  The  alpha  angle  describes  the  predominant  alignment  of  the  debris  in 
the  channel  with  respect  to  the  How  direction  These  indices  were  first  used  by  ('berry  and 
Beschta  (1086)  in  connection  with  I,WD  flume  experiments.  The  Beta  angle  is  a  measure  of 
the  predominant  orientation  of  the  debris  jam  material  in  the  vertical  plane,  at  ninety  degrees 
to  the  flow  direction.  See  figure  3. 1  below 

Figure  3.1  Alpha  and  Beta  Indices  (modified  from  Cherry  and  Beschta,  1986) 

Alpha  Index  (degrees)  Beta  Indes  (degrees) 


5)  Sinuosity  :  A  visual  estimate  of  channel  sinuosity  in  the  jam  reach  (straight,  slightly  sinuous, 
sinuous,  meandering). 

b)  Knickzoncs  :  Presence  of  Kickpoints  or  knickzoncs  in  the  channel  reach  (a  measure  of 
channel  instability) 

7)  Sediment  :  Bedload  I)  ,,  classification  (clay,  silt,  sand,  gravel) 

8)  Deposition/Scour  :  Estimated  volume  of  bar  deposition,  backwater  sedimentation  and 
bed/bank  scour  induced  by  each  jam.  Total  deposition  and  scour  values  are  then  calculated  for 
each  channel  reach 

The  debris  volume  and  debris  frequency  measures,  morphological  classification  jam 
sedimentation  and  erosion,  and  net  jam  sediment  bubgets  in  each  reach  have  been  and  plotted 
against  three  independant  catchment  variables.  These  are;  reach  drainage  basin  area,  average 
reach  channel  width;  and  reach  average  unit  stream  power.  Unit  stream  power  is  calculated 
using  the  following  equation: 


CO 


(3.1) 


pgQs 
w 

where:  (o  stieam  power  per  unit  bed  area  N/m/s  • 

p  density  of  water  kg/W 
g  gravitational  constant  (0.81  m/s/s) 

Q  predicted  Q2  discharge  (cumecs)  q 

s  =  bed  slope  (m/m) 
w  -  reach  average  channel  width  (m) 

These  independent  catchment  variables  are  used  to  determine  whether  the  geomorphological 

4p> 

effects  of  LWD  have  a  coherent  and  predictable  spatial  relationship 
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Figure  3.2  Geological  map  of  Mississippi  showing  the  DEC  survey  area  (modified  from 
Watson  et  al.,  1993) 


Figure  3.3  DEC'  project  area  site  location  map  (mollified  from  Raphelt  ot  al,  1905) 


4  RESULTS 

Results  are  divided  into  analysis  of  debris  input  Mechanisms.  Processes,  and  Debris  jam 
residence  times 

4.1  INPUT  MECHANISMS 

Figure  4. 1  shows  a  table  of  debris  input  mechanisms  at  each  jam  site  for  each  reach  which 
have  significant  debris  jams.  The  dominant  input  mechanism  is  through  outer  bank  erosion  in 
active  meanders  (43%),  followed  by  input  due  to  reach  scale  instability  (30%).  Random  input 
processes  (windthiow,  beaver  activity  and  floated  material)  accounted  for  only  37%  of  debris 
input  in  total.  Debris  input  in  these  unstable  rivers  is  therefore,  in  general,  caused  by  spatially 
predictable  phenomena  Debris  input  zones  and  major  dam  accumulation  can  thus  be  located, 
without  resorting  to  expensive  and  time  consuming  field  reconnaissance,  by  using  secondary 
map  data  to  locate  channel  reaches  which  have  wooded  riparian  zones  with  active  meandering 
and/or  lateral  instability  due  to  degradation.  In  catchments  up  to  about  50  miles  square  jams 
appear  to  form  where  the  key  debris  (that  is  large  whole  trees  which  initiate  jam  formation)  fall 
into  the  channel.  Jams  are  therefore  commonly  located  in  bend  apices  or  in  unstable  reaches 
downstream  of  knickpoints.  Figure  4.2  demonstrates  this  observation,  showing  debris  jam 
locations  just  downstream  of  bend  apices  on  a  planform  plot  of  Abiaca  Creek.  Jams  do  not, 
however,  appear  to  have  a  regular  or  predictable  spatial  distribution  relating  to  channel  bed 
topography,  such  as  a  location  on  riffle  heads 

4.2  IMPACTS  UPON  CHANNEL  MORPHOLOGY 

4.2.1  Impact  upon  bed  topography 

Comparison  of  thalweg  plots  from  reaches  which  are  located  in  wooded  riparian  /.ones  with 
those  whicn  arc  located  in  agricultural  riparian  zones  show  that  bed  topography  is  generally 
much  more  varied  in  wooded  reaches.  The  difference  in  morphology  must  be  due,  in  the  most 
pail,  to  the  presence  of  debris-induced  bed  scour  and  sedimentation  because  the  other 
controlling  variables,  including  sediment  load  and  substrate  geology  are  identical 
This  difference  can  be  illustrated  by  comparing  thalweg  plots  for  Redbanks  ('reek  (figure 
4.40),  which  has  a  puiely  agricultural  riparian  zone,  with  Worsham  (’reek.  West  Fork  (figure 
4.3 1 )  which  has  a  wooded  riparian  zone 

4.2.2  Debris  Jam  Reconnaissance  analysis 

Refer  to  figure  4  3  showing  the  geoinorphological  reconnaissance  results  foi  the  May  1005 
survey. 
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The  following  geomorphological  variables  have  analysed  and  plotted  against  three  independent 
catchment  variables;  drainage  basin  area,  stream  power  and  channel  top  width;  to  determine 
whethei  the  geomorphological  effects  of  LWD  have  a  coherent  and  predictable  spatial 
relationship. 

1)  Number  of  Jams  per  unit  reach  length 

Figure  4.4  shows  a  plot  of  the  number  of  debris  jams  per  1  OOOtf  of  channel,  for  each  reach 
surveyed,  against  drainage  basin  area.  A  negative  relationship  was  expected  in  thi.  plot 
because  of  the  greater  flow  competence  in  channels  with  a  larger  catchment,  however,  no  r'ear 
fun  ion  is  evident  in  the  plot  Jam  frequency  has  also  been  plotted  against  reach  average 
composite  channel  top  width,  and  reach  unit  stream  power.  These  are  graphed  in  figures  4.5 
and  4.6.  Again,  however,  no  relationship,  positive  or  negative  is  evident  in  either  plot.  A 
simple  spatial  relationship  is  therefore  too  simple  an  approach  to  understanding  debris  jam 
dynamics  and  it  must  be  recognised  that  debris  jam  frequency  (and  therefore  volume)  in 
particular  reaches  is  also  heavily  dependant  on  the  comparative  rate  of  debris  input  processes 
at  work,  in  balance  with  transport  from  the  reach.  Input  rates  may  vary  enormously,  especially 
in  the  environment  in  question,  where  heavy  debris  input  through  bank  erosion  is  dependant 
upon  channel  instability  and  active  meandering  processes,  factors  which  are,  to  a  certain 
degree,  i  ndependent  of  the  reach  catchment  area. 

?)  Volume  of  debris  per  unit  reach  length 

Figure  4.7  shows  a  plot  of  the  volume  of  debris  present  in  each  channel  per  1 00011  equivalent 
ag  inst  drainage  area.  A  negative  trend  in  this  plot  might  have  been  expected  because  more- 
debris  is  likely  to  be  mobilised  in  larger  catchments  due  to  greater  flow  competence  and  the 
increasing  channel  dimensions  as  compared  to  the  average  key  debris  size.  There  is  no 
significant  trend  evident,  however,  the  plot  being  almost  horizontal,  with  the  exception  of  one 
or  two  reaches,  most  notably  Abiaca  6  (drainage  area  00  mi;)  It  is  possible  that,  rather  than 
debris  volume  being  a  smooth  function  of  drainage  area,  a  step  in  the  debris  transport  capacity 
exists  due  to  some  threshold  discharge  value  for  transport  of  the  key  in-ehannel  debris. 
Although  a  step  exists  between  the  volume  in  Abiaca  0  anil  the  other  reaches  in  smaller 
catchments,  there  is  insufficient  data  ovei  the  full  range  of  drainage  area  displayed  to 
demonstrate  such  a  phenomena  conclusively.  Plots  ol  debris  volume  against  composite  channel 
width  (figure  4  H)  and  stream  power  per  unit  channel  widlh  (figure  T  ■>)  also  show  no 
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significant  directional  trend,  again  where  negative  relationship  may  have  reasonably  been 
expected. 

3}  Flow  deflection  classification 

Figure  4.10  shows  a  bar  chart  of  jam  morphological  classification  plotted  against  drainage 
basin  area  classes  with  a  10  mi'’  interval.  Unfortunately,  because  there  are  not  an  equal  number 
of  i  ‘hservations  in  each  category  it  is  very  difficult  to  make  a  meaningful  interpretation  of  the 
different  distributions  in  each  category.  If,  however,  the  categories  containing  debris  jams  are 
taken  as  representational  of  all  the  jams  present  over  that  range,  the  distribution  shows  a  weak 
trend  from  dam  type  jams  in  catchments  less  than  10  mi7  through  to  a  dominance  of  parallel 
jams  in  larger  catchments  (>  40  mi7). 

4)  Flo  vy  blockage  classification 

The  influence  that  each  jam  has,  in  terms  of  blocking  flow  and  sediment  movement  is 
qualitatively  described  using  the  Active/Complete/  Partial  classification  (see  chapter  3).  Partial 
jams  make  up  82%  of  the  total,  complete  jams  14%  and  active  jams  only  4%  (see  figure  4.3). 
The  majority  of  jams  are  therefore  partial  and  do  not  cause  significant  backwater  effects  and 
backwater  sediment  wedge  formation.  Log  steps,  which  are  a  majoi  feature  associated  with 
debris  jams  in  gravel  bed  rivers  (see  Macdonald,  1982)  are  therefore  not  a  prevalent  feature  in 
these  unstable,  sand-bed  creeks. 

5)  Alpha  angle 

The  Alpha  variable  describes  the  mean  angle  of  orientation  of  debris  relative  to  the  channel 
flow  direction.  There  is  no  distinct  relationship  with  increasing  drainage  basin  area  (no  plot 
shown),  although  it  could  be  theorised  that  the  alpha  angle  would  increase  from  90  to  180 
degrees  as  drainage  area  (and  therefore  channel  width  and  mean  discharge)  becomes  larger 
due  to  greater  How  forces  in  larger  channels  causing  the  rotation  of  trees  parallel  with  the 
current.  The  results  in  figuie  4.3  show  that  filly  eight  percent  of  the  sites  have  a  mean  debris 
direction  of  90  degrees,  24  percent  have  a  mean  direction  of  180  degrees,  and  18  percent  have 
a  mean  direction  between  90  and  180  degrees.  It  is  evident  therefore  that  the  majority  of  trees 
enter  the  channel  at  90  degrees  to  the  flow,  some  arc  subsequently  rotated  by  the  fo'ce  of 
flow,  but  these  observations  demonstrate  that  the  tree  lengtli/channel  width  index,  which  has 
been  used  to  describe  blockage,  and  flow  disturbance,  by  debris,  is  quite  justified  for  modelling 
purposes  (see  Wallenstein  1995). 

9)  Met  a  angle 


V) 


The  Beta  variable  describes  the  mean  angle  of  debris  in  relation  to  the  channel  bed.  Debris  that 
rests  on  the  channel  bed  along  its  entire  length  has  a  value  of  0  degrees,  while  debris  which  has 
one  end  in  the  channel  and  the  other  end  supported  on  the  bank  has  a  Beta  value  ranging 
between  0  and  90  degrees.  Referring  to  figure  4.3,  data  show  that  eighty  five  percent  of 
debris  jams  have  an  average  beta  angle  of  0  degrees,  while  only  15  percent  have  a  beta  angle 
greater  than  5  degrees.  Thus,  at  the  majority  of  sites  debris  will  offer  maximum  blockage  area 
at  low  flows,  helping  to  retain  sediment  transported  in  the  flow,  but  at  the  same  time  causing 
flow  constriction  and  deflection  which  may  result  in  bank  failure.  Debris  inclined  more  than  5 
degrees  appears  to  be  more  prevalent  at  creeks  which  have  a  low  width  to  depth  ratio,  such  as 
Nolehoe  creek  where  the  w/d  ratio  approaches  1  in  certain  sections,  the  obvious  reason  for 
this  relationship  being  that  trees  falling  into  these  channels  will  rotate  more  than  90  degree 
before  coming  to  rest  (unless  they  span  the  channel). 

Debris  Beta  angle  will  affect  the  morphological  impact  that  a  jam  has  upon  the  channel. 
Steeply  inclined  debris  is  unlikely  to  cause  significant  flow  deflection  and  will  have  a  low 
sediment  retention  capacity  as  compared  with  debris  lying  parallel  to  the  channel  bed. 
Significant  bank  erosion  or  bar  formation  is  therefore  unlikely,  however,  localised  scour 
around  the  debris  will  cause  pools  to  form  akin  to  those  caused  by  flow  round  bridge  piers. 
Debris  lying  on  the  channel  bed  along  its  full  length  will  have  greater  interaction  with  the  flow 
and  a  range  of  features  may  develop  depending  upon  the  sediment  type,  and  stage  of  flow. 
Such  debris  inclination  also  offers  at  greater  upstream  area  on  which  to  trap  floating  debris, 
and  these  jams  will  consequently  grow,  and  become  more  coherent  structures,  at  a  faster  rate 
than  inclined  debris  formations.  See  figure  4  !  I 


Figure  4, 1 1  :  Flow  patterns  around  debris  with  different  inclinations 
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7)  Sedimentation  and  scour  associated  whh  debris  jams 

Sedimentation  and  scour  values  tor  each  site  are  shown  in  figure  4.3.  These  values  have  been 
graphed  in  figure  4.12  against  drainage  basin  area,  against  unit  channel  width  in  figure  4  13, 
and  against  stream  power  per  unit  width  in  figure  4  14  Triangles  represent  sediment  retention, 
through  backwater  and  bar  deposits,  for  each  creek,  averaged  over  1000  ft.  Rectangles 
represent  bed  scour  due  to  flow  restriction  and  turbulence  around  debris  jams  for  each  creek, 
averaged  over  1000  tt  of  channel  While  no  distinct  pattern  emerges  in  any  of  the  three  plots, 
in  figure  4  12  it  is  evident  that  sedimentation  and  scour  values  are  low  at  sites  less  than  10  mi7, 
then  increase  to  a  peak  at  between  20  and  40  mi7  and  then  appear  to  decline  in  larger 
catchments.  This  trend  is  also  vaguely  apparent  in  figure  4.13  with  sedimentation  rising  to  a 
peak  at  a  channel  width  of  around  58  tf  ihen  declining  to  zero  at  108  ft  The  peak  value  for 
scour  in  figure  4  12,  discounting  the  extreme  value  at  24  mi7  ,  lies  at  around  10  mi7  while  the 
peak  value  for  sedimentation  lies  at  around  42  mi2  This  peak  distribution  is  conti adicted 
somewhat  in  the  composite  width  plot  where  scour  values  peak  at  a  larger  width  than  that  for 
sedimentation,  however,  the  general  pattern  can  be  explained  if  it  is  linked  to  the  idealised 
deb'  i  channel  interaction  curves  outlined  in  Wallerstein  (1904'  k  values  for  bed  scour  lie 
in  the  zone  where  underflow  and  dam-type  jams  are  likely  to  be  prevalent,  that  is,  where  flow 
is  constricted  and,  in  the  underflow  dam  case,  forced  to  impinge  on  the  channel  bed  resulting 
in  the  formation  of  scour  holes  Scour  values  then  tail  away  as  flow  constriction  around  debris 
is  reduced  as  channel  widths  become  greater  than  the  average  debris  length  Sedimentation 
values  peak  at  slightly  larger  drainage  areas,  where  jams  arc  either  of  the  dam  or  deflector 
type  Dam  type  jams  will  trap  sediment  the  backwater  area,  while  deflector  dams  often  have 
a  bar  downstream  of  the  jam  due  to  energy  dissipation  and  low/mid  flow  deflection  towards 
the  channel  bank  In  larger  channels  still,  jams  become  flow  parallel  red  occupy  a  smaller 
percentage  of  the  channel  cross  section,  so  that  sedimentation,  through  energy  dissipation 
around  the  debus,  been  ...s  less  •elective  The  idealised  distribution  of  debris  related 
sedimentation  and  scour  with  increasing  drama}*.-*  basin  area  1:;  shown  in  figure  4.15  below  it 
must  be  recognised  t hat  it  ;s  often  difficult  to  attribute  an  is  of  sedimentation  or  scour 
specifically  to  a  debris  jam  T’  is  is  especially  the  ease  with  areas  of  bar  sedimentation  in  larger 
channels,  where  a  ciicular  explanation  pioWctu  arises  in  i ho i  liars  may  either  he  initial  d  by 
dingy  dissipation  ">  the  Ice  of  a  debris  build- up,  or  may  pre-date  the  debiis,  which 
subsequently  comes  to  icvi  upon  its  upsticam  face  this  problem  wre;  recognised  by  llickin 
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(1084)  when  trying  to  explain  th.  formation  of  scroll-bars  In  either  case,  however,  debris 
would,  at  the  very  least,  enhance  sediment  deposition  on  the  bar 

Figure  4.15  :  Idealised  diagram  of  debris  related  sedimentation  and  scour 
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The  magnitude  of  the  scour  and  sedimentation  peaks  may  well  vary  with  bedload  and 
suspended  sediment  type,  with  coarse,  gravel,  loads  causing  the  sedimentation  peak  to  rise, 
and  the  bed  scour  peak  to  fall  proportionately,  and  vice-versa  with  a  fine  sand  or  silt  sediment 
load 

The  LWD  sediment  budget  for  each  reach  surveyed  (averaged  over  1000  ft  of  channel),  that 
is,  the  difference  between  the  total  sedimentation  volume  and  total  scour  volume,  induced  by 
debris,  for  each  channel  reach  is  plotted  against  drainage  basin  area  in  figure  4.16  against 
composite  channel  width  in  figure  4.17  and  against  stream  power  per  unit  channel  width  in 
figuie  4.18  Relation  to  die  three  independent  catch"'  "t  w.rkf  'ler  reveals  no  recognisable 
spatial  trends  However,  quite  importantly,  the  graphs  do  show  that  the  balance  between 
sediment  scour,  and  sediment  deposition  caused  by  debris  jams  is  in  favour  of  net 
sedimentation,  with  nine  of  the  channel  reaches  having  a  positive  budget  and  only  loui  having 
a  negative  budget  The  total  difference  between  positive  and  negative  budgets,  summed 
together  for  all  channel  reaches  leaves  a  net  positive  balance  of  68  m'  of  sediment.  This  value 
would  be  much  higher  if  it  were  not  tor  the  one  anomalous  negative  budget  value  of -260  m\ 
This  value  relates  to  Harland  Creek  (reach  ! )  and  is  thought  to  be  partly  due  to  overestimation 
of  debris  related  scour  in  channel  bends  where  it  was  hard  to  distinguish  between  meander 
pools  and  scour  induced  by  flow  constriction  around  debris.  The  net  positive  sediment  balance 
lends  weight  to  the  argument  that  debris  lias  an  overall  beneficial  impact,  in  terms  of  sediment 
retention  in  these  highly  unstable  channels  with  excessive  sediment  yields,  although  the  total 


retention  values  are  possibly  negligibly  small  when  compared  with  total  sediment  yields  from 
each  reach 

4.3  DEBRIS  JAM  RESIDENCE  TIMES 

Figure  4. 19  to  4.40  show  thalweg  plots  for  each  reach  with  surveyed  positions  of  each  debris 
jam  site  marked  on.  Data  from  the  May  1995  survey  are  overlaid  onto  the  May  1994  plots  so 
that  an  assessment  can  be  made  of  the  rate  of  debris  input  and  of  the  longevity  of  jams,  which 
determines  their  effectiveness  as  geomorphological  channel  controls. 

Initially,  it  was  hoped  that  consecutive  years  of  bed  survey  data  could  be  overlaid  to  monitor 
changes  in  topography  around  each  jam  site.  Unft.  rtunately,  however,  it  has  been  found  that 
there  are  too  many  data  survey  errors  to  allow  a  realistic  comparison  of  local  scale  features 
between  the  two  years  of  data. 

Thalweg  plots  are  not  available  for  Coila  Creek,  Abiaca  Creek  (site  3)  and  Hic'kahala  Creek 
(site  22)  due  to  severe  errors  in  the  data. 

The  following  creeks  have  debris  jam  sites  surveyed  into  the  long  profile  data  : 

Worsham  C reek  (Wr  .l  Fork)  Harland  Creek  ( 1 )  Hickahala  Creek  (11) 

Worsham  Creek  (Middle  Fork)  Lick  Creek  Nolehoe  Creek 

Sykes  C reek  Fannegusha  Creek  Abiaca  C reek  (6) 

Lee  Creek  Long  Creek  Perry  Creek 

Of  these  the  following  reaches  has  jam  sites  surveyed  in  both  May  1994  and  May  1995 
Harland  Creek  Hickahala  Creek  Lick  Creek 

Nolehoe  Creek  Abiaca  Creek  (4) 

The  above  plots  give  an  indication  of the  stability  of  jams  and  the  rate  of  debris  input  over  the 
one  year  peiiod 

Harland  Creek  ( I)  :  Figure  4  19 

Two  major  jams  (I  and  2)  surveyed  in  1994  were  still  intact  in  1995.  In  addition,  six  new  jam 
sites  wore  located,  al1  of  which  were  caused  by  tree  topple  due  to  bank  instability 
Hickahala  Creel  Figure  4  2G 

bite  I  had  been  removed  by  1995  due  to  the  construction  of  a  new  road  bridge.  Sue  2,  a  jam 
formed  from  branches  l  ro'iyii!  down  in  an  ice  storm,  had  been  washed  down  to  jam  3  by  may 
1995  Site  4  remained  intact  in  1995  and  had  collected  a  large  quantity  of  fresh  material. 


Lick  Creek  Figure  4.21 

Site  1  remained  intact  between  1994  an  t  1995.  Five  new  underflow  type  jams  with  trunks 
resting  across  the  channel  at  bankfull  level  were  surveyed  in  1995.  Many  of  these  are  unlikely 
to  interfere  with  the  flow  except  in  the  most  extreme  flood  events.  Sites  2,  3  and  4  had  been 
removed  by  May  1995  due  to  the  constn  :ion  of  a  new  high  drop  grade  control  structure. 
Nolehoe  Creek  :  Figure  4. 22 

All  nine  sites  surveyed  in  1994  were  found  to  be  intact  in  1995.  Several  had  gained  new 
material  both  floated  down  from  upstream  and  from  tree  topple  due  to  flow  deflection  and 
bank  undercutting  around  the  existing  jams.  This  reach  is  highly  unstable  which  accounts  for 
the  very  large  volume  of  debris  present. 

Abiaca  Creek  (site  no.  4) :  Figure  4  23 

Sites  1  and  2,  surveyed  in  1994  were  found  to  be  intact  in  1995  These  jams  have  been  formed 
through  key  debris  input  where  the  outer  bank  of  two  active  meanders  have  eroded 

It  is  evident  from  the  comparison  of  this  one  year  interval  data  that  the  majority  of  jams  have 
remained  intact  over  this  period,  while  a  number  of  new  jams  have  formed  New  jams  appear 
to  have  been  initialed  either  though  bank  failure  in  unstable  reaches,  such  as  on  Harland  Creek 
or  through  a  combination  tree  blow  down  on  unstable,  over  steepend  banks,  such  as  on  the 
upper  reaches  of  l  ick  Creek 

C'onclusions  about  short-term  jam  stability  can  only  be  made  when  the  survey  data  is 
compared  with  coi  responding  discharge  data  for  ihc  one  year  interval.  This  will  show  the 
magnitude  and  frequency  of  events  that  the  jams  are  able  to  survive.  Unfortunately,  as  yet  the 
relevant  data  has  not  been  made  available  to  permit  this  analysis.  Given  the  warm,  humid 
climate,  and  highly  credible  nature  of  these  channels  it  is  likely,  however,  that  jams  arc  more 
transient  features  in  northern  Mississippi  than  in  Pacific  North  West  gravel  bed  rivers  because 
debris  will  decompose  quickly  and  become  bypassed  and  transported  downstream  more 
readily.  A  long-term  monitoring  programme  is  required,  however,  to  verify  this  asseifion 

New  debris  jam  sites  were  surveyed  at  six  creek  reaches  It  is  not  known  if  these  jams  were 
present  in  May  1994  because  time  constraints  meant  that  debris  reconnaissance  could  not  be 
carried  out  at  these  reaches. 


Sykes  Creek  :  Figure  4  24 

Four  sites  were  surveyed  into  this  reach  Key  debris  input  is  primarily  the  result  of  outer  bank 
failure  in  active  meanders. 

Fannegusha  Creek  :  Figure  4. 25 

This  reach  is  highly  unstable  and  has  just  had  a  grade  control  constructed  at  it’s  downstream 
end  Three  major  debris  jams  were  located  within  the  reach,  all  cause  by  bank  failure  input. 
Abiaca  Creek  (site  6)  :  Figure  4.26 

This  reach  is  situated  at  the  point  where  the  creek  flows  out  of  the  Bluff  Line  hills  onto  the 
Mississippi  floodplain.  The  gradient  is  shallow  and  the  channel  is  agrading  slightly.  The  two 
jams  present  appeared  to  have  been  formed  by  material  that  has  floated  down  from  upstream 
rather  than  from  local  tree  input. 

Lee  Creek  :  Figure  4.27 

Four  small  jams  were  found  in  this  reach.  Two  in  the  upstream  section  have  been  created  by 
debris  input  due  to  bank  instability,  while  the  two  lower  down  the  reach  appear  to  have  been 
formed  by  blown  down  trees 
Long  Creek  :  Figure  4.28 

This  entire  reach  is  choked  with  debris,  and  “palaeodebris”,  that  is  woody  material  that  must 
have  been  preserved  and  buried  in  the  alluvium  of  older  channels,  and  is  now  being  exhumed 
by  the  creek  as  it  erodes  into  it’s  bed  and  banks  Distinct  jams  are  very  hard  to  define  in  this 
reach,  but  nine  concentrations  of  debris  were  surveyed  into  the  thalweg  data.  This  reach  was 
highly  unstable  in  the  past  but  is  now  beginning  to  stabilise  due  to  the  construction  of  a  major 
grade  control  at  the  downstream  end  of  the  reach,  and  teveral  smaller  structures,  upstream 
through  the  reach  section 
Pen  v  Creek  :  Figure  4.2() 

This  r  each  displays  all  the  phases  of  the  channel  evolution  sequence  (see  Schutnm,  Harvey  anti 
Watson,  I '>84),  from  type  5  at  the  lower  limit  of  the  site  to  undisturbed  type  I  at  the  upstream 
end  Seven  jam  sites  were  located  in  this  reach  Two  jams  in  the  lower  reach  have  been  created 
by  blown  down  tra  three  in  the  upper  third  of  the  reach  are  located  in  unstable,  type  2, 
channel  ant!  have  it  .'died  from  bank  instability,  and  ihc  final  two  jams,  located  in  type  1 
channel,  in  the  upstream  end  of  the  reach,  have  been  created  by  bend  outer  bank  erosion. 


The  following  two  reaches  had  jams  surveyed  in  May  1994  but  were  not  survey  in  May  1995 
owing  to  time  constraints: 

Worsham  Creek  (Middle  Fork) :  Figure  4.30 

Thirteen  sites  were  surveyed  into  this  reach.  Debris  input  had  occurred  mainly  through  bank 
failure,  and  the  entire  reach  appeared  to  be  degrading  and  laterally  unstable. 

Worsham  Creek  (West  Fork)  Figure  4.5 1 

Seven  sites  were  surveyed  into  this  reach.  Sites  four  and  five  were  beaver  dams,  while  the 
other  jams  were  created  by  debris  input  through,  slab  type,  bank  failures. 

The  remaining  reaches  (figures  4.32  to  4.40)  have  no  debris  jams  surveyed  into  the  thalweg 
plots.  Otoucalofa  Creek  (figure  4.33),  James  Wolf  Creek  (figure  4.34),  Sarter  Creek  (figure 
4.36),  Abiaca  Creek,  site  21  (figure  4.37),  Burney  Branch  (figure  4.38),  Harland  Creek,  site 
23  (figure  4.39)  and  Red  Banks  Creek  (figure  4.40)  were  all  surveyed  for  debris  but  were 
found  to  contain  no  jams.  The  reason  for  an  absence  of  debris  material  was  found  to  be  due  to 
three  factors.  Either,  one,  because  the  reach  is  stable  and  has  no  major  bank  erosion  (Sarter 
Creek,  Burney  Branch  and  Abiaca  (’reek,  site  21),  two,  because  the  reach  has  no  woody 
riparian  zone  to  supply  material  to  the  channel  (James  Wolf  Creek  and  Red  Banks  Creek),  or, 
three,  because  the  discharge,  and  channel  dimensions  of  the  creek  are  such  that  even  large 
whole  trees  which  fall  into  the  creek  do  not  remain  in  one  location  long  enough  to  establish  a 
coherent  jam  before  being  swept  away  downstream  (Otoucalofa  Creek,  Harland  Creek,  site 
23,  and  Abiaca  ('reek,  site  2 1 ). 

Worsham  ('reck.  East  Fork  (figure  4.32)  and  llotopha  ('reek  (figure  4.35)  have  not  been 
suiveyed  for  debris  jams,  but  it  is  suspected  that  debris  is  present  in  these  reaches  owing  to  the 
fact  that  both  are  known  to  be  degrading  and  may  well  therefore  have  unstable,  failing,  banks 


Figure  4.1  LWD  Input  Mechanisms  (May  1995  Survey) 
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ABIACA  CREEK  (SITE  4) 

Planform  Survey  :  May  1994  /  1995 


anform  oiot  of  Abiaca  Creek.  Site  No. 


Figure  4.3  LWD  Geomcrphological  Reconnaissance  Data  {May  1996)  £ 
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Jam  Frequency  against  Drainage  Area 

Reach  average  values  :  May  1995 
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Figure  &£  Jam  frequency  against  drainage  basin  area  plot 


Jam  Frequency  against  Channel  Width 

Reach  average  values  :  May  1995 
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Figure  4.5  Jam  frequency  against  composite  channel  width  plot 


Jam  Frequency  against  Stream  Power 

Reach  average  values  :  May  1995 
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“igu-e  4.6  Jam  frequency  against  unit  stream  power  piot 


Debris  Volume  against  Drainage  Area 

Reach  average  values  :  May  1995 
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i-  gjre  4.7  Deans  volume  against  drainage  basin  area  plot 


Debris  Volume  against  Channel  Width 

Reach  average  values  :  May  1995 
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Composite  Channel  Width  (ft) 


Debris  Volume  against  Stream  Power 

Reach  average  value  :  May  1995 


aaainst  unit  stream  cower  plot 


Jam  classification  /  drainage  area 

May  1995  data 


Drainage  Basin  Area  Ranges  (MiA2) 


Debris  Induced  Sedimentation  &  Erosion 

Reach  average  values  :  May  1995 
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debris  related  sedimentation  ■  debris  related  bed  scour 


Debris  Induced  Sedimentation  &  Erosion 

Reach  average  values  :  May  1995 


~i_-e  ^.'3  Debris  induced  sedimentation  and  erosion  :  composite  width  plot 


Debris  Induced  Sedimentation  &  Erosion 

Reach  average  values  :  May  1995 
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debris  induced  sedimentation  and  erosion  :  unit  stream  power  plot 


VD  Sediment  Budget 

average  values  :  May  ^99 


sedimentation  -  scour 


LWD  sediment  budget 

Reach  average  values  :  May  1995 
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LWD  sediment  budget :  composite  channel  width  plot 


LWD  sediment  budget 

Reach  average  values  :  May  1995 
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HARLAND  CREEK 

Thalweg  Survey  Profile 


pro'l'e  Dlot  of  Harland  Creek.  Site  No.  1  :  May  1994/95  data. 


LICK  CREEK 

Thalweg  Survey  Profile 


o 

o 

o 

LO 


o 

o 

o 


o 

o 

o 

00 

-4— « 

<D 

<D 


o 

o 

o 

CM 


4 


^  <o 

t 


o 

o 

o 


I 


o 

o 

o 

o 

CM 

T — 

o 

CD 

CO 

CO 

CO 

CM 

}00J 


o 


6 


hi 


CD 

~i 

CD 

LL 


■  Debris  Jam  Site,  May  1995  a  Debris  Jam  Location,  May  1994 


NOLEHOE  CREEK 

Thalweg  Survey  Profile 
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ABIACA  CREEK  (SITE  4) 

THALWEG  SURVEY  PROFILE 
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■  Debris  Jain  location.  May  1995  a  Debris  Jam  location,  May  1994 


SYKES  CREEK 

Thalweg  Survey  Profile 
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B  Debris  Jam  location  May  1995 


FANMEGUSHA  CREEK 

Thalweg  Survey  Profile 
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re  4.25  Long  profile  plot  of  Fannegusha  Creek  :  May  1995  data 


ABIACA  CREEK  (SITE  6) 

Thalweg  Survey  Profile 
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■  Debris  Jam  location.  May  1995 


LEE  CREEK 

Thalweg  Survey  Profile 


■0  4.2?  ^ong  prove  plot  of  i_e9  Creek  :  May  1995  data 


LONG  CREEK 

Thalweg  Survey  Profile 


|  Debris  Jam  location.  May  1995 


PERRY  CREEK 

Thalweg  Survey  Profile 
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■  Debris  Jam  location.  Mav  1 995 


WORSHAM  CREEK  (MIDDLE  FORK) 

Thalweg  Survey  Profile 
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Debris  Jam  location,  May  1994  data 


WORSHAM  CREEK  (WEST  FORK) 

Thalweg  Survey  Profile 


Debris  Jam  location.  May  1994 


OTOUCALOFA  CREEK 

Thalweg  Survey  Profile 
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gure  4.33  Long  profile  plot  of  Otoucalofa  Creek  :  May  1994  data 


HOTOPHA  CREEK 

Thalweg  Survey  Profile 


-cure  4.35  L:'C  profile  plot  of  Hotopfa  Creek  :  May  1994  data 


SARTER  CREEK 

Thalweg  Survey  Profile 
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Igure  4.35  Lcng  profile  plot  of  Sarter  Creek  :  May  1994  data 


ABIACA  CREEK  (SITE  21) 

Thalweg  Survey  Profile 


ig-'e  4.37  Long  profile  plot  of  Abiaca  Creek,  Site  21  :  May  1994  data 


BURNEY  BRANCH 

Thalweg  Survey  Profile 
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Figure  4.38  Long  profile  plot  of  Burney  Branch  :  May  1994  data 
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Figure  4.39  U 


5  COMPUTER  MODELS 

5.1  LARGE  WOODY  DEBRIS  MANAGEMENT  PROGRAM 

The  relationships  between  LWD  formations  and  channel  processes  have  been  incorporated 
into  a  LWD  Management  Program.  Version  1 .0  of  this  program  is  included  with  this  report  on 
a  disk.  Appendix  1  a  user  manual  for  this  program. 

This  is  an  updated  version  of  the  program  to  that  which  was  included  with  the  Project  R&D 
7258-EN-09,  submitted  to  the  US  Army  Corps  of  Engineers,  June  1995  (Wallerstein,  1995). 
This  program  predicts  the  jam  type  in  a  given  reach,  determines  its  impact  upon  the  channel 
and  outlines  an  appropriate  management  strategy.  Inputs  variables  are  those  which  have  been 
found  to  be  most  critical  in  the  LWD  system  and  include  channel  width  (determined  from  a 
catchment  area  function),  average  riparian  tree  height,  reach  sediment  type  and  the  riparian 
landuse  type.  The  ratio  of  tree  height  to  channel  width  is  used  to  define  the  debris  jam  type 
present,  with  the  precise  limits  of  each  classification  determined  from  the  empirical 
relationships.  Sediment  diameter  is  used  to  give  an  indication  of  the  jams  potential  to  induce 
backwater  sedimentation  or  downstream  bars.  Debris  jam  types  are  classified  using  a  scheme 
modified  from  Robinson  &  Beschta  (1990),  described  in  Wallerstein  &  Thome  (1994).  Jam 
types  are  divided  into  Underflow,  Dam,  Deflector  and  Flow  Parallel.  Figure  5.1  shows  this 
classification  scheme.  The  program  output  takes  the  form  of  a  text  file  which  describes  the 
classification  chosen,  and  offers  basic  in-channel  LWD  management  strategies.  While  the 
management  strategies  are  based  solely  on  theoretical  considerations,  the  program  never-the- 
less  provides  a  framework  for  future  model  development  as  empirical  relationships  between 
the  variables  are  better  characterised.  A  flow  diagram  of  the  computer  program  is  shown  in 
figure  5.2. 

The  program  has  also  been  linked  to  a  GIS  (Geographical  Information  System)  front  end 
which  has  been  constructed  by  Peter  Cheeseman,  a  masters  student  at  Nottingham  University 
(see  Cheeseman,  1995).  The  project  was  carried  out  to  demonstrate  the  potential  for  using 
GIS  as  a  platform  for  data  input  to  expert  systems  to  aid  engineers  with  river  basin 
management. 

The  GIS  was  constructed  in  ARC  Info  using  data  layers,  supplied  by  the  WES  Intergraph 
data-base,  for  the  Abiaca  Creek  watershed  and  provides  automatic  data  input,  for  the 
necessary  variables  and  a  platform  for  running  the  program.  This  watershed  was  selected 
because  it  contains  four  debris  survey  reaches  which  are  being  monitored  in  the  current 
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research.  The  theoretical  model  can  therefore  be  tested  against  the  empirical  data  results  from 
the  field  studies,  and  be  validated  and  further  developed.  This  management  model  is  simple, 
but  provides  a  framework  for  future  development  as  empirical  relationships  between  variables 
are  better  characterised. 


Figure  5.1  Revised  Debris  Jain  Classification  Scheme 


Figure  5.2  LWD  management  program  flow  diagram 
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5.2  DKHKIS  AT  BRIDCK  PIER  PREDICTION  PROGRAM 

This  piogram  calculates  the  probability  of  debris  build-up  at  bridge  piers,  ami  the  associated 
debris  induced  scour,  based  upon  modified  theoretical  equations  published  by  Melville  and 
Dongol  (1992)  and  Simons  and  l.i  (1979).  Version  1.0  of  this  program  is  also  included  on  the 
enclosed  disk.  Appendix  2  contains  a  user  manual  foi  this  piogram 

5.2.1  Bridge  Scour  Model  with  Debris  Build-up 

There  are  only  a  limited  number  of  studies  that  have  addressed  the  problem  of  debris 
accumulations  at  bridges  Melville  &  Dongol  (1992)  look  at  the  problem  of  pier  scour  due  to 
debris,  while  .unions  &  Li  (see  Callander  1980)  have  used  a  probabilistic  approach  to  quantify 
the  rate  of  bridge  span  blockage  by  debris  and  the  subsequent  backwater  effect  and  pressure 
forces  generated  on  the  piei  . 

Local  scour  at  bridge  piers  h  is  been  extensively  investigated.  However  the  impact  of  debris 
rafls  at  piers  which  create  additional  llow  obstruction  and  therefore  increase  scour  depths  has 
been  largely  neglected.  A  design  method  for  estimation  of  scour  depths  at  piers  is  presented 
by  Melville  and  Sutherland  (1988),  based  on  envelope  curves  from  laboratory  data.  The  hugest 
local  scour  depth  at  a  cylindrical  pier  is  estimated  to  he  2  41)  where  I)  is  the  pier  diameter. 
This  vaiue  is  reduced,  however,  using  multiplying  factors  where  dear-water  scour  conditions 
exist,  the  flow  is  relatively  shallow,  and  lire  sediment  size  relatively  coarse.  In  the  ease  of  lion- 
cylindrical  piers,  additional  multiplying  factors  are  apph  d  to  account  for  piers  shape  amt 
alignment  Consideration  of  the  likelihood  of  debris  build-up  is  not  addressed  by  Melville  and 
Dongol  ( 1992)  but  they  do  note,  however,  that  single  cylindrical  piers  are  the  least  likely  shape 
to  accumulate  debris,  and  that  tire  free  space  betv  :cn  columns  is  soMom  great  enough  1o  pass 
debr  is.  Prediction  of  the  size  of  possible  debr  is  rail  accumulations  remain  ;  i  h*  biggest  problem 
for  accurate  factor  of  safety  calculations. 

The  experimental  arrangement  used  by  Melville  and  Dongol  is  shown  in  figure  5..L 
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Figure  S.3  Experimental  Set-up  (modified  T  o  in  Melville  and  Dongol,  1992). 


The  design  curve  tor  pier  scour  without  debris  accumulations,  developed  by  Melville  and 
Sutherland  is  described  by  the  following  two  equations: 
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—  1.872 
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(S.D 


f/v 
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-2.4 


>2.6 


(5.2) 


This  shows  that  scour  depth  increases  with  increasing  llow  depth  towards  a  limiting  value  for 
Y/P  -2.6.  The  same  trend  is  found  for  piers  with  debris  accumulations  for  valtK  s  of  Y/I>  4  At 
higher  values  of  Y/l)  scour  depths  decrease  again  because  the  proportion  of  pier  length 
covered  by  debris  decreases  For  deep  Hows  the  effect  of  debris  would  become  insignificant 
and  tend  towards  If.  value  ds/l)  2  4 


The  effective  mametci  of  a  pier  with  a  debris  accumulation,  l)e,  is  given  by, 

'V  '  (5  i) 

y 

According  to  equation  T  1  l>c  is  calculated  as  a  weighted  average  of  an  c  libel  ivc  length  'I'd 
0  ST'I'd  of  i he  debris  rail  with  diametci  Pd  and  a  length  of  the  piei  (Y  I'd  )  vvitli  diameter  I) 
(see  figure  5.1).  flic  factoi  0.12  was  dcln  mined  hv  evaluating,  the  linn ■  s  of  Td  and  Pd/P  for 
the  hypothetical  case  when'  P  is  assumed  to  he  /cm  md  the  dehris  is  :<  ..aimed  In  extend  to  the 


base  of  the  scour  hole 


0 


SO 


I)  can  therefore  be  substituted  for  !)e  to  calculate  scour  depth  at.  piers  with  debris 
accumulations  using  the  Melville  and  Sutherland  design  method  Conversely  a  maximum 
allowable  Td  and  Dd  can  be  calculated  by  specifying  an  upper  scour  depth  within  an 
acceptable  factor  of  safety  for  a  given  pier  size. 

5.2.2  Probability  based  debris  build-up  model 

The  rate  of  debris  accumulation  at  a  bridge  is  difficult  to  quantify.  The  only  method  found  in 
the  literature  is  that  presented  by  Simons  Li  (107b)  in  an  MSc  thesis  by  Callander  (1980) 
entitled  "fluvial  processes  occurring  at  bridge  sites”. 

According  to  Simons  £  Li,  the  trapping  efficiency  of  a  bridge  is  determined  by: 

1 )  Clearance  beneath  the  bridge 
Span  lengths 

.! )  Size  and  concentration  of  debris  elements 

flic  following  possible  consequences  are  identified  which  can  result  from  debris  blockage: 

1 )  Backwater  effects 

2)  Potential  local  How  diversion 

3)  Channel  avulsion 

4)  Bridge  failure 

Simons  &  Li  express  the  volume  of  dehris  as  a  fraction  of  the  sediment  yield,  and  suggest  a 
vegetation  debris  yield  of  1%.  In  an  attempt  to  estimate  the  number  and  volume  of  trees 
arriving  at  a  bridge  they  utilise  the  volume  of  flood-plain  erosion  necessary  fo  yield  a  tree,  and 
use  a  representative  free  size  for  the  watershed. 

'Lees  are  assumed  u>  be  cylindrical  with  a  diameter  !)t,  and  a  height  fit.  The  span  between 
piers  is  L  >  and  the  clearance  between  the  water  sm  lace  and  the  under  side  of  the  bridge  is  C. 
The  chance  that  a  tree  will  he  (nipped  depends  on  a  larger  diameter  however,  l)b,  which 
represents  either'  the  canopy  dimension  or  i'»e  root  /.oik  .  whichever  is  larger  See  figure  5.4. 

Ill  If  Ls  I  lie  probability  of  at  least  one  average  tree  being  trapped  is  100%.  The  blocked  area 
is  then  estimated  to  be.  NlltHl,  where  N  is  the  equivalent  number  of  average  trees  assumed  to 
be  trapped  against  the  upstream  face  of  the  bridge 
lfl  It  Ls  a  probabilistic  approach  is  used 

PI  is  the  probability  of  a  tree  being  trapped,  and  as  the  blockage  beneath  a  span  increases  so 
the  chance  of  ollici  trees  being  happed  inc-ieases  The  probability  of  the  first  tree  being 


trapped  is  assumed  to  b  i  ratio  of  half  the  tree  diameter,  Db,  to  the  total  waterway  area 
beneath  a  span,  LsC 

Oh1  /A)  7T  Db' 


PTl  =- 


(5.4) 


(5.0) 


LsC  8  LsC 

(,i  (see  Callander,  1080)  observed  that  a  tree  caught  on  a  pier  will  in  general  lie  with  its  trunk 
in  the  direction  of  flow.  A  tree  thus  trapped  offers  an  area  of 

X(rr/t/>2/4)=rr/8/J/>;  (5.5) 

to  trap  other  debris. 

In  general  when  (m-l)  trees  arc  trapped  beneath  a  span  the  probability  of  an  mth  tree 
becoming  trapped  is 

!  8 

LsC  — (m  —I  )(7r/ ■>/>'’  /  8) 

The  probability  of  passing  all  NT  trees  from  the  watershed  is 
( 1-PT!  )NI 

The  probability  of  at  least  one  tree  being  trapped  at  a  span  is 
PI  =  I-(I-PT1)n 

where  N  is  the  equivalent  number  of  average  trees  arriving  at  the  span.  According  to  l.i  most 
trees  will  stay  close  to  the  bank,  thus; 

N  NT/2  (5.9) 

The  probability  that  m  trees  will  be  trapped  is; 

Pm  [  I  -( I  -PTm)|N‘,I"’lllJP(ni- 1 )  (5.10) 

On  this  basis  the  probability  of  a  least  m  trees  being  trapped  (for  any  m  ■  N)  can  be  estimated. 
In  order  to  calculate  'I’d  and  Dd  there  needs  to  be  an  estimate  of  the  blockage  area  It  is 
assumed  that  debris  elements  stack  up  and  that  trees  overlap  by  Ob/2  Thus  for  m  trees 
trapped  the  percentage  of  the  waterway  area  which  is  blocked  is 

/h(],k  Dir  j A) 


(5.7) 

(5.8) 


%li  lockage  = 


Ls( ' 


xi  oo% 


(5.11) 


Having  er  >  1  ed  m  and  knowing  l)b  (he  increase  depth  of  water  (wd)  at  the  bridge  is 
assumed  to  be 

Air</  J ml )//  /  2  (5  12) 

The  blockage  generates  a  prcssuie  force  (PI)  which  ae.ts  normal  to  ti  e  bridge  is 

/’/■  fry. ml /  4  O  I  H 


1 


(s  1  t) 


where  y  is  the  specific  weight  of  water 

Figure  5.4  Schematic  diagram  of  debi  uniulation  at  bridge  piers 


5.2.3  Program  construction  and  input  varial  s 

In  essence  the  program  runs  the  Simons  and  Li  probability  model  and  then  calculates  tire 
potential  pier  scour  due  to  a  debris  mat  size  based  upon  the  blockage  area  assuming  all  the 
trees  available  in  the  reach  upstream  become  trapped  beneath  a  span. 

Initially  average  (mid  height)  tree  trunk  diameter  (I)t),  a  maximum  tree  diameter  (l)h),  (eii 
root  wad  or  canopy,  whichever  is  the  larger),  and  average  tree  height  (111)  values  are  entered 
Next,  the  number  of  trees  approaching  the  bridge  span  (NT)  is  entered  Although  Simons  and 
l.i  suggest  using  N  NT/  hi  the  probability  calculations,  this  model  assumes  that  all  the  trees 
available  in  the  upstream  reach  will  pass  through  the  span  in  question  I  lowevei .  the  number  of 
spans  (S)  between  piers  (P)  that  are  set  in  the  channel  will  normally  he  S  PM  (counting  the 
two  spans  between  pier  and  river  hank)  It  is  therefore  ncccssaiy,  for  an  accurate  assessment 
ot  blockage  potential  and  debris  related  scour,  to  calculate  probabilities  for  each  span 
individually,  pci  Imps  using  a  simple  division  rule  (N  NT  /  S)  for  N  trees  ariiving  at  each 
span.  It  is  let)  up  to  tlx-  user  to  make  the  appropriate  adjustments  for  each  span  NT  can  either 
he  estimated  in  the  field  ami  entered  as  a  total  potential  tree  supply  or  can  he  estimated 
thiougli  calculation  of  potential  hank  failuit  in  the  upslicsm  reach  To  calculate  the  latter 
estimate  a  lipaiian  lice  density  value  is  icquiied.  the  length  ol  the  icucli  in  question  and  the 
potential  bank  1  hline  width  Mr  liiilnie  width  value  can  he  detei  mined  using  an  appmpiiale 


bank  stability  model  such  as  BURBANK  (Burgi,  1095).  The  potential  number  ol'trees  that  will 
reach  the  span  is  then  calculated  as; 

tree  density  x  failure  width  x  reach  length  x  2  (two  banks) 

Finally,  the  bridge  pier  diameter  (D),  span  between  piers  (Is)  and  average  (low  depth  (Y) 
values  are  entered 

Calculations  then  proceed  as  follows: 

1 )  If  tree  height  is  less  than  the  pier  spacing  the  probability  of  the  first  tree  becoming  caught  is 
calculated,  followed  by  the  probability  of  the  next  tree  becoming  caught  consecutively  This  is 
repeated  for  n  trees  up  to  NT. 

In  the  calculation  of  trapping  potential  it  is  considered  that,  the  use  of  the  ratio  of  tree  area  to 
the  entire  area  under  the  span  as  suggested  by  Simons  and  Li,  is  somewhat  inappropriate  as 
tree  capture  is  dependant  only  upon  the  length  of  span  and  diameter  of  tree  given  that  the 
water  level  is  constant  Deck  elevation  above  the  water  ((')  has  therefore  been  substituted  with 
maximum  tree  diameter  (Db)  in  this  model. 

2)  If  tree  height  is  greater  than  span  width  it  is  assumed,  as  outlined  in  the  Simons  and  Li 
model  that  at  least  one  tree  will  become  trapped  and  thus  all  subsequent  trees  arriving  at  the 
span  will  also  be  caught. 

3)  The  percentage  of  the  channel  cross  sectional  area  that  is  blocked  if  all  the  trees  supplied  to 
the  reach  become  trapped  is  calculated  as  outlined  in  the  theoretical  model  (see  figure  5.4)  for 
I  It  •  Is.  1  lowever  if  I  It  -  Is  !)t  is  substituted  for  Db  and  the  blockage  area  is  calculated  as  : 
(((square  root  x  blockage  area)  blockage  depth  (assuming  debris  builds  up  as  a  square)  x 
tree  height)  /  (span  width  x  (low  depth))  x  100  % 

This  calculation  assumes  that  for  lit  •  Is  all  trees  will  build  up  in  a  square  formation,  but  at  00 
degrees  to  the  llow  direction,  as  oppose  to  parallel  with  the  How  when  I  It  •  Is. 

4)  The  hydrostatic  pressure  force  on  each  pier  per  unit  width  is  calculated  as  : 

pressure  force  bulk  weight  of  water  x  blockage  depth  x  I  (unit  width)  x  (blockage  depth  / 

2). 

5)  Itiidgc  pita  scorn  with  the  debris  accumulation  is  then  calculated  using  the  Melville  and 
Dongnl  model.  If  lit  ■  Is  the  debris  rail  diameter  is  taken  as  the  square  root  of  the  blockage 
area  (assuming  debris  build-up  is  in  a  square).  If  III  •  Is  debris  iafl  diameter  is  assumed  to  be 
lit  because  the  debris  is  aligned  parallel  with  the  direction  of  How  The  scour  depth  is 
calculated  using  the  base  value  of  2.41)  If  the  additional  multiplying  factor  :,  arc  required  they 


H4 


must  be  added  by  hand,  referring  to  the  graphs  and  tables  supplied  in  Melville  and  Dongols’ 
1988  paper.  Factors  which  reduce  (his  value  are  applied  where  elear-water  scour  conditions 
exist  (K|),  the  flow  depth  is  relatively  shallow  (Kv),  and  the  sediment  is  relatively  large  (K,i).  If 
pieis  are  not  cylindrical  two  additional  factors  arc  required;  a  shape  factor  (K.,)  and  an 
alignment  factor  (K„).  The  following  data  are  required  to  calculate  these  additional  factors  : 
mean  approach  velocity  for  the  design  flood  (I  J);  median  particle  size  (d5o);  standard  deviation 
of  the  particle  distribution  (o„)  dm/dw;  pier  diameter  (D),  angle  of  flow  attack;  pier 
dimensions;  and  pier  shape, 

It  should  be  noted  that  the  formula,  developed  by  Melville  and  Dongol  for  calculating  debris 
related  pier  scour,  was  only  developed  for  floating  debris  accumulations.  However,  it  is 
considered  by  the  authors  that  this  formula  can  be  extended  to  debris  accumulations  which 
have  their  base  resting  on  the  channel  bed,  as  (lie  critical  factor  in  the  calculation  method  is  an 
effective  pier  diameter,  which  is,  in  any  ease,  extended  to  the  channel  bed  in  the  situation 
where  the  debris  is  floating. 


6  CONCLUSIONS 


The  following  conclusions  can  he  drawn  from  the  results  obtained  in  this  study  : 

1)  The  sources  of  major  debris  input  to  the  channel  network  can  be  predicted  by  mapping  the 

areal  distribution  of  reaches  with  wooded  riparian  zones  which  are  also  laterally 
unstable  or  actively  migrating. 

2)  The  impact  of  debris  jams  in  sand-bed  rivers  is  different  to  that  in  gravel  bed  rivers  as 

distinct  log  steps  do  not  form. 

3)  Debris  is  a  key  factor  controlling  channel  bed  topography  in  sand-bed  rivers,  creating  a 

more  heterogeneous  profile  than  found  in  debris  free  reaches.  This  is  important  for 
creating  productive  aquatic  habitat. 

4)  Debris  jams  in  sand-bed  rivers  are  stable  in  the  short  term  (I  yr .).  But  total  residence  times 

may  be  longer  in  gravel-bed  rivers. 

5)  Reconnaissance  evidence  suggests  that  jam  forms  change  in  a  predictable  manner 

downstream  through  the  channel  network. 

b)  Debris  jams  dissipate  flow  energy,  and  reduce  sediment  routing  rates,  and  therefore  do  not, 
exacerbate  channel  bed  degradation  problems.  The  disliibution  of  sedimentation  and 
scour  .  ssoeiated  with  debris  jams  appeals  have  an  "comorphically  explainable 
distribution  when  related  to  drainage  basin  area 

mana<;kiwknt  rkcommknd  tions 

The  folk)'.  •!>:■  management  iceoinmeiulation:;  liavi  been  made  based  upon  the  findings 
discuss'  >1  herein: 

1)  Basin- wide  debt  is  clearance  is  linn-  ccssary  ami  may  have  doliimental  al  loots  D.  In  i:;  does 

not  :mpcai  to  .aceilsitc  the  degradation  pioblem  encountered  in  the  bluff-line  sti earns 
and  can  even  accelerate  sggiadstion  processes  locally. 

2)  Debris  should  be  loll  in  place  or  even  introduced  into  the  channel  if  the  management  aim  is 

habitat  enliaiu  ement 

I)  Debus  jams  may  have  In  he  removed  ovci  the  mid  range  ol  basin  sizes  if  the  piimaiy 
management  aim  is  to  i educe  hank  ciosion 

4)  In  rivers  with  urn  of-nvci  stiiiclures,  excessive  debus  load  >,  may  be  paitially  udueed 
through  eoiiiml  ofiipslic.am  clianiK  I  dcgiadalioa  and  oiius  bank  emsiou 
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APPKNDIX  I 

Large  Woody  Debris  Management  Program  (Version  1)  liser  Manual 

Introduction 

This  manual  describes  version  I  of  an  on-going  developmental  program  designed  to  aid 
engineers,  geomorphologists  and  planners  with  management  of  woody  debris  in  river  channels 
throughout  the  catchment  network. 

The  executable  program  for  this  model  (DEBRIS)  is  written  in  C++  and  is  contained  on  the 
supplied  MS-DOS  formatted  disk,  labelled  “DEBRIS  &  DPROB”  The  file  called  debris.cpp 
contains  the  program  source  code 

This  program  may  be  freely  copied,  but  the  source  code  may  not  be  distributed  to  other 
parties.  The  author  accepts  no  responsibility  or  liability  resulting  from  the  use  of  the  program. 

Installation 

To  install  DFBRIS  and  DPROB  insert  the  supplied  disk  into  the  A  drive  and  type  a  install.  A 
directory  called  LWD  will  then  he  created  on  the  C  drive  and  the  files  debris  exe  and 
dprob  exe  copied  into  it 

Running  the  program 

It  should  be  noted  that  the  required  units  for  data  entry  in  this  program  are  SI  (millimetres, 
metres,  kilometres). 

Run  the  program  by  typing  “debris”  from  the  I.WD  directory. 

1 )  F.nter  a  file  name  for  the  output  text  tile  e  g  “basin  1  txf” 

2)  Enter  the  riparian  vegetation  type  that  is  predominantly  found  along  the  reach  of  river  in 
question.  A  simolc  distinction  is  made  between  wooded  and  agricultural.  Enter  “w”  or  “a” 

If  “agricultural  is  entered  the  next  prompt  will  ask  if  you  wish  to  run  the  program  again  Type 
“y"  or  “n” 

4)  If  “w”  is  entered  the  next  prompt  will  ask  you  to  enter  the  average  riparian  tree  height 
(metres). 

5)  Entei  ihe  drainage  basin  area  (kilometres  square)  of  the  reach  in  question  This 
parameter  is  used  to  calculate  an  average  channel  width  based  upon  a  selected  width  function. 
<i)  Average  channel  width  can  either  be  calculated  using  the  default  function,  derived  by 
Schumm,  Harvey  and  Watson  (l()84),  which  is  appropriate  for  channels  in  northern 
Mississippi,  or  the  user  can  enter  their  own  function  The  default  function  is  . 

average  channel  width  |(4(>.77  x  I  >a" ' ')  t  ( 1 5.7  x  Da",w  )|  /  2 


'll 


Where  Da  -  drainage  basin  area  in  square  miles  This  is  converted  to  square  kilometres  by  the 
program 

Select  “n”  to  use  the  default  function  or  select  “y”  to  enter  a  new  function  A  new  function 
takes  the  form  of 

average  channel  width  -  constant  (a)  x  drainage  area  cona,*n,  ,h) 

The  user  is  prompted  first  for  constant  a,  then  for  constant  b  both  of  which  must  he 
numerical  values  from  functions  developed  in  St  units  Note  constant  b  must  be  less  the  I  else 
the  program  will  crash 
7)  Enter  the  channel  bed  sediment  (mm) 

9)  Tite  user  will  then  be  asked  whether  they  wish  to  run  the  program  again.  Enter  “y”  or  “n”. 
If  “y”  is  yped  the  program  will  mn  again  and  return  to  the  prompt  asking  for  the  riparian 
vegetation  type.  If  “n”  is  selected  the  program  will  terminate  and  return  to  the  operating 
system. 

Program  Output 

The  results  are  written  to  an  output  file  which  has  the  file  name  specified  by  the  user  The 
output  file  is  created  in  the  LWD  directory  If  the  program  is  unable  to  create  the  output  file 
an  error  message  will  appear  reading  '  Error  opening  file”.  If  this  happens  the  program  should 
be  terminated 

The  output  text  file  can  be  viewed,  once  the  program  has  been  ended,  in  any  text  editor  or  in  a 
word  processor  such  as  Word  or  Word  Perfect 

Successive  outputs  from  each  program  will  be  added  sequentially  to  the  same  text  file. 

Program  runs  are  distinguished  by  a  numbered  header  (i.e.  program  run  number  1 . 

program  run  number  2 .  etc  ) 

Results  have  the  following  standard  outpu*  format: 

1 )  A  program  run  and  “  Woody  Debris  Management  Output”  header 

2)  A  list  of  the  input  parameters  including  the  calculated  channel  width. 

3)  A  description  of  the  debris  jam  classification  type  chosen  and  the  geomorphologica!  impact 
of  debris  jams  for  that  reach.  Management  recommendations  (displayed  in  capital  letters) 
based  upon  the  geomorphological  processes 

l  est  C  ases 

These  examples  arc  provided  to  <usuie  that  the  user  is  familiar  with  the  procedures  used  to 
mn  the  program,  and  to  check  t hat  the  supplied  program  is  working  correctly 


1 .  Start  the  program  by  typing  “debris”  from  the  l.WD  directory 

2)  Enter  an  output  fde  name  e  g.  “basin l.txt” 

3)  Enter  the  following  values  : 


run  1 

run  2 

run  3 

run  4 

riparian  vegetation 

w 

w 

w 

a 

average  tree  height 

10 

10 

10 

- 

drainage  basin  area 

40 

200 

20 

- 

alternative  width  function  9 

n 

n 

y 

- 

constant  a 

- 

- 

10 

- 

constant  b 

- 

- 

0.30 

- 

sediment  type 

0.25 

1.0 

2.6 

- 

7)  Leave  the  program  and  run  a  text  editor  or  word  processor  to  view  the  output  file. 

8)  T  e  output  for  runs  1  to  4  should  be  as  follows; 


Program  run  number  I 

W(X)DY  Ui'BRIS  MANAtil-MliNT  OUTPUT 


ITe  riparian  land  lype  in  the  chosen  reach  is  wixxlcd 
Hie  drainage  basin  area  is  40  kiloinctics  square 
Idle  average  ripai'iun  tree  height  is  1(1  metres 
Ilie  sediment  I  )50  in  this  reach  is  0.25  min 

C  hannel  width  calculated  using  the  Schinnm,  llaivev  &  Watson  formula 
The  average  channel  width  in  this  reach  is  20.7052  metros 

(lid  >Mt)RIT  It  )I ,( X JIC'AI ,  IMPACT  OK  WvX)|)Y  DlllRIS  IN  I  t  IP  Slil  .l  XTId)  KKAC'I  I 


Average  debris  length  is  less  tham  ().(>  times  the  cluimiel  width 

Large  Wixxly  I  >ehris  ( I.WI ))  input  mechanisms  to  this  reach  will  include 

floatation  of  debris  from  upstream  reaches,  tiee  toppple  due  to  bunk  litiluie 

m  unstable  reaches,  and  also  due  to  hank  failure  in  b  aid  apices  if  the 

channel  in  this  reach  is  meandering  Although  failed  trees  arc  likely  to  cntc-i 

the  channel  at  W  degrees  to  the  channel,  Hows  will  lx;  compctant  enough  to  rotate 

debris  so  that  MX >W  PARAI.I.1.1 . I  .WD  debris  puns  form.  I.WI )  will  also  lx*  lrans|X)ited 

downstream  is  high  llows,  and  deposited  against  the  bank-base  on  the  outside  of 

meander  bends  or  at  bridge  pieis  and  other  mn-of-the-rivci  structures. 

(Sec  Wnllerstem  &  Thome  ( l'W4),  for  a  complete  description  of  1 IX W-PAKAI.I ,1(1.  tvjy.i  jams) 
lb-  negative  geomotphological  impact  of  this  lyjK-  of  lam  ill  terms  of  bank 
ci  'ion  and  bed  scour  is  minimal  Hank  toes  may  even  lx  stabilised  by  debris  build-up 
I XX iris  mnv  also  initiate  or  aeelerate  the  fonnation  ol  inid-clianni  I  and  lateral  bars 

WOODY  l)KI IRIS  MANAOl-MTNT  KP.Ct  'MMI  NI  )A  I  It  )NS 
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DbBRIS  Cbb.ARANOt  I INNIiCbSSARY  II-  II'  IS  KI-.YI  IN To  PI  ACli  AT  HANK  IOI  S  OR  IN  BARS. 

URI1X.1S  ANROTHHR  RIJN-OT-TI  Ili-RIVliR  STRUC TUKliS  Si  IOI  II  OBI!  MONITORltD  I  OR  JAM  BIJII.D-UP 
(See  Melville  A  Dongo!  ( I  992)  for  a  method  of 'calculating  piei  scorn  due  to  debris  build-up 
and  Callandar  ( 19X0)  for  a  probabilily  based  tuelhiHl  of  predielng  debris  build  up  al  bridge 
piers,  and  formula  for  calculating  lire  increased  pressure  force  and  How  allies  due  lo  debris) 

moating  dihris  may  aiso  provi:  to  hr  a  dangi-r  tosmai.i  boats 


Rill  HRKNCbS 

CAbbANUliR  ( 1980),  fluvial  processes  oecuring  ai  bridge  sites, 

MSc.  Thesis,  Colorado  State  University. 

MlibVlbl.li  &  IX)NCK)1,  ( 1992),  Bridge  Pier  Scour  with  1  )ebris  Accumulation, 

Journal  of  Hydraulic  Hngineering,  ASCIi,  vol  I  IX,  no.  9 

SCI  HI  MM.  1 1AR  VKY  A  WATSON  ( 19X4),  Incised  Channels  Morpliologv.  Dynamics  and  Control, 
Water  Resources  Publication,  pp,  1 1 1-1 S9 

WAbbliRSTlilN  A  l  l  IORNH  ( 1 994)  Impact  of  in-ehaimel  organic  debris  on  fluvial  process 
and  channel  morpliologv,  Yazoo  Basin,  Mississippi,  t  Iniversitv  of  Nottingham,  Department 
of  Geography,  working  pa|ier  no.  29 


Program  mil  number  2 

W(X)1)Y  I  >I1BR1S  MANAObMI'NT  <  )l  I  I  PI  I  I 


'I Tie  riparian  land  type  in  the  chosen  reach  is  wooded 
The  drainage  basin  area  is  200  kilometres  square 
The  average  riparian  tree  height  is  10  metres 
'Hie  sediment  1)50  in  tlii.v  reach  is  I  inni 

Channel  width  calculated  using  the  Sctunnm.  I  laivev  A  Watson  fonnula 
Ihc  average  channel  width  in  this  reach  is  5Y2b2<>  incites 

edit  )M(  )RPI  It  )I  ,(X  HCAI  IMPACT  Of  W(X)|)Y  |)l  BRIS  IN  1 1  Hi  SliblXTKD  RliA<  II 

Average  debris  length  is  less  tliam  ().(>  times  (lie  channel  width 

barge  Woody  I  XT  :  (l,WI ))  input  inech  'iiisms  lo  tbu  reach  will  include 

lloataliou  of  debus  (nun  upstream  teaches,  lice  loppple  due  to  hank  lailuic 

in  unstable  readies,  and  also  due  to  bank  failure  in  bend  apices  if  Hie 

channel  ill  this  reach  is  meandering.  Although  tailed  trees  me  likely  lo  eiitei 

the  channel  at  90  degrees  to  ihc  channel.  Hows  will  tie  competant  enough  to  rotate 

debris  so  that  Id ,()W  PARAI  I  lil.  I ,W1 )  debris  jams  lonu.  I  AVI)  will  also  Ire trans)iorted 

downstream  is  high  Hows,  and  deposited  against  the  bank  base  on  the  outside  of 

meander  bonds  or  at  bridge  piers  and  other  nm-of-the  liver  stnictures. 

(See  Wallerstein  A  'Phonic  ( 1994),  for  a  complete  description  ol  fbt  1W-PARA1  .bill .  type  jams) 

I  lie  negative  geonioipliological  impact  of  this  tvjie  of  jam  in  terms  of  ‘sink 
erosion  and  bcrl  scour  is  minimal  Hank  toes  mav  even  lie  stabilised  bv  debris  build-up 
I  >cbns  iiuiv  also  initiate  or  aceleiate  Hie  fra  nialion  of  inid-charmel  and  l  iteral  bars. 

W(K)I)Y  HUBRIS  MANAGIiMliNT  KliCOMMb.NI  )A  1  It  )NS 

I  >|-.BRIS  (  I  .IIARANCb.  1  INNI  CbSSARY  If  I  f  IS  Kb  Yl-.l )  IN  I  < )  PI  ACI-  A  I  HANK  TOPS  OK  IN  BARS. 

BRIIXilS  AND  Oil  IliR  RUN  ( I  I'-'l  I  Ib-RIVIR  S  IRIK  II  IRIS  SI  |OI  l|.|  1  HI  I  MONITt  IRIiD  H)R  JAM  BIJII.D  UP 
(See  Melville  A  I  krngol  ( 1 9‘>? )  for  a  mclluxl  of  calculiiliiii'  piei  .coni  due  lo  debt  i::  build-up 
and  t  allaiidar  ( 1 9X0)  lor  a  probability  based  method  of  pn-di-tu;;  debris  build  up  at  budge 


'Mi 


piers,  and  formula  for  calculating  the  inci eased  pressure  force  and  How  afflux  due  tu  debris) 
H.OATINU  DidiRIS  MAY  AI.SO  PROVh  d  O  lib  A  DANOLR  l'( )  SMALL  1IOATS 
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I’rognun  nut  tuunhet  3 

W(X  )l )Y  DLURIS  MANAOLMHNI  OlllPIJ'l 


lire  riparian  land  type  in  the  chosen  reach  is  wooded 

Hie  drainage  basin  area  is  20  kilometres  squire 

The  average  riparian  tree  height  is  10  metres 

Hie  sediment  DM)  in  tins  reach  is  2.6  min 

Channel  width  calculated  using  the  following  formula 

width  =  (drainage area  '  0.3)  x  10 

Hie  average  channel  width  in  this  reach  is  24  5(>4<>  metres 

(ILOMORITIOl.tKdCAI.  IMPACT  Oh  W(X)I  >Y  DI  BRIS  IN  I  I lb  SLI.LC  I'Ll)  RLACI I 

Average  debris  length  is  less  tham  0  <i  limes  the  channel  width. 

large  Wendy  Debris  (l.VVD)  input  mechanisms  to  this  reach  will  include 

lloalation  of  debris  fiom  upstream  renehes,  bee  loppple  due  to  bank  failure 

in  unstable  reaches,  and  also  due  to  brink  failure  in  bend  apices  if  the 

channel  m  this  reach  is  mcandcrir"  Although  failed  trees  ate  likely  to  entei 

the  channel  at  40  degrees  to  the  channel .  Hows  will  tic  compctmil  enough  to  totalc 

debt  is  so  that  il. I  )W  I’AKAI  d.LL  I  ,WI )  debris  jams  form  I.WI )  will  also  be  tiansporicd 

downstream  is  high  tlows.  and  deposited  against  the  hank  base  on  the  outside  of 

meander  bends  or  at  budge  piers  and  other  mii-ol-'lie-iiver  slliictmes. 

(Sec  Wullerstein  it  ihorne  ( 1444 1,  for  a  complete  description  of  M.( )W-1’ARAI  I.LL  t\|ie  jams) 

Hie  negative  geomorphologieal  impact  of  this  type  of  gun  in  tennx  of  hank 

erosion  and  bed  scour  is  minimal  Hank  Iocs  may  even  lie  stabilised  by  debris  build-up 

Debris  may  also  initiate  oi  ncelerute  the  Ibnnation  of  mid  channel  and  lateral  bars 

WtXlDY  DI  BRIS  MANAOIMI.N  I  KbCOMMI  NDA  LIONS 

DLliRIS  CI.LARANCI.  I INNLCLSSAKY  If  I  I  IS  KLYId )  IN  I (>  PI  ACL  Ad  BANK  IOLS  ( )|<  IN  BARS 
IIRIIXILS  AND  Odd II  R  RUN  Od-I'llb -RIVLK  S  I'RUC'I  ! IRIiS  SI  It  MILD  HL  MONI  lOKLI d  I  OR  JAM  BUILD  III 
(See  Melville  it  1  dongol  ( 1 442)  lor  a  luclhod  of  calculating  pier  scout  due  to  debus  build-up 
and  Callandai  (14X0)  for  a  probability  based  method  ofprcdii  mg  debris  build  opal  bridge 
piets,  and  for  inula  for  calculating  the  increased  presume  (iirce  and  (low  afflux  due  to  delu  is) 

I  l.(  )A  I  INO  DldlRIS  MAY  Al  SO  I’ROVL  TO  lib  A  HANOI  K  TO  SMAI.I.  lit  )ATS 


CALLANDER  (1980),  Eluvial  processes  occnring  at  bridge  sites, 

MSc.  'Diesis,  Colorado  State  University, 

MLLVILLE  &- 1 X )N(iOI  (1992),  Bridge  Her  Scour  witli  Debris  Accumulation, 

Journal  ot  Hydraulic  Engineering,  ASCH,  vol.  1 '8,  no  l) 

SCHUMM,  HARVKY  <fe  WATSON  (19)14),  Lncised  Channels  :  Morphology,  Dynamics  and  Centro1 
Water  Resources  Publication,  pp.  111-159 

WALLERSTK1N  &  THORNE  (1994),  Impact  of  in-dmnne)  organic  debris  on  lluvial  process 
and  channel  morphology,  Yazoo  Basin,  Mississippi,  University  of  Nottingham,  Department 
of  Ueographv,  working  paper  no  29 


Program  run  number  4 


VVtXTDY  DEBRIS  MANAGEMENT  OUTPUT 


The  riparian  land  tyjie  in  the  chosen  retieh  is  agrieultur.il 


Because  the  inunediate  riparian  /one  is  agricultural  land,  ot  the  reach 
selected  is  in  open  water  woody  debris  input  through  uec  blowdown  or 
death,  beaver  activity  and  channel  bank  erosion  will  he  minimal 
Substantial  debris  jams  are  therefore  unlikely  to  build  up  in  the  retieh 

W(X)DY  DHBRIS  MANAGEMENT  RECOMMENDATIONS 

DEBRIS  JAMS  NOT  PRESENT  IN  REACH,  DHBRIS  REMOVAL  UNKCHSSARY 

sm.vli.hr  debris  may  be  input  from  upstream  reaches  however  and  debris 

BUILD-UP  AT  STRUCTURED  SUCH  AS  BRIDGE  PIERS.  UX'KS,  DAMS  AND  WIERS,  SHOULD  HE  MONITORED. 

IF  CHANNEL  BED  IS  UNSTABLE  AND  HAS  A  SAN  I). SEDIMENT  LOAD  CONSIDERATION 

SHOULD  BE  GIVEN  TO  ARTIFICIAL  DEBRIS  INPUT  ON  THE  EHQUENCY  OF  THE  EXPECTED 

RIFFLE  SPACING  TO  VARY  CHANNEL  BED  TOPOGRAPH  Y  AND  T  HUS  IMPROVE  AQUATIC  1 IAHITAT 

AND  EN1 JANCE  CHANNEL  STABILISATION  THROUGH  SEDIMENT  RETENTION 

LARGE  W(X)DY  DHBRJS  CAN  ALSO  HE  KEYED  INTO  THE  OUTSIDE  ERODING  BANK  OF  ACTIVE 

MEANDER  1)1  iNIXS  TO  REDUCE  LOCAL  NEAR-BANK  SHE  AR-ST  RESS 
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APPENDIX  2 

Debris  at  Bridge  Pier  Prediction  Program  (Version  1.0)  User  Manual 

introduction 

This  manual  describes  a  probability  based  model  for  predicting  the  build-up  of  debris  at  bridge 
piers  It  is  based  upon  two  theoretical  models  developed  by  Melville  and  Sutherland  (1992) 
and  Simons  and  Li  (in  Callander,  1980). 

The  executable  program  for  this  model  (DPROB)  is  written  in  C++  and  is  contained  on  a  MS- 
DOS  formatted  disk,  labelled  “DEBRIS  &  DPROB”.  The  file  called  dprob  cpp  contains  the 
program  source  code. 

This  program  may  be  freely  copied,  but  the  source  code  may  not  be  distributed  to  other 
parties.  The  author  accepts  no  responsibility  or  liability  resulting  from  the  use  of  the  program. 

Installation 

To  install  DEBRIS  and  DPROB  insert  the  supplied  disk  into  the  A  drive  and  type  a  install  A 
directory  called  LWD  will  then  be  created  on  the  C  drive  and  the  files  debris.exe  and 
dprob  exe  copied  into  it 

Running  the  program 

The  required  units  for  data  entry  in  this  program  are  SI 
Run  the  program  by  tying  “dprob"  from  the  LWD  directory. 

1)  Enter  a  Tile  name  for  the  output  text  file  e  g  “results  1  txt”. 

2)  Enter  the  average  riparian  tree  trunk  diameter  (m) 

3)  Enter  the  average  maximum  root  wad  or  canopy  diameter  (whichever  is  larger),  (m). 

4)  Enter  the  average  riparian  tree  height  (m) 

5)  The  number  of  trees  approaching  the  bridge  span  can  be  either  be  entered  as  a  total 
value  from  upstream  reach  observations,  or  can  be  estimated  through  the  input  of  a  riparian 
tree  den.'ty,  bank  top  failure  width  and  reach  length  Type  “I  for  an  estimated  value.  Type 
“2"  to  enter  an  observed  value 

5a)  If  “2”  is  typed  enter  the  number  of  trees  available  in  the  upstream  reach. 

5b)  If  “I”  is  typed 

i)  Enter  average  riparian  tree  density  (tn/m) 

ii)  Enter  bank  top  failure  width  (m)  This  can  be  calculated  using  an  appropriate  bank 
stability  model 

iii)  Enter  reach  length  (m) 


b)  Enter  bridge  pier  diameter  (m) 

7)  Enter  distance  between  bridge  piers  at  water  level  (mV 

8)  Enter  flow  depth  immediately  upstream  of  bridge  (m). 

9)  Enter  lly”  to  input  new  values  or  enter  “n”  to  end  the  program. 

Program  Output 

The  results  from  DPROB  are  written  an  output  file  which  has  the  file  name  specified  by  the 
user  during  the  first  program  run  The  output  files  will  be  created  in  the  LWD  directory  If  the 
program  is  unable  to  create  the  output  file  an  error  message  will  appear  reading  “Error 
opening  file”.  If  this  message  appears  the  program  should  be  terminated. 

The  output  text  files  can  be  viewed,  once  the  program  is  terminated,  in  a  text  editor  or  word 
processor  such  as  Word  or  Word  Perfect. 

Successive  outputs  from  each  program  run  will  be  added  sequentially  to  the  same  text  file 
Program  runs  are  distinguished  by  a  numbered  header  (i.e.  progiam  run  number  1;  program 
run  number  2; . etc  ) 

Results  have  the  following  standard  output  formal 

1 )  A  “Debris  at  Bridge  Pier  Calculations”  header. 

2)  A  list  of  the  input  parameter  values,  including  the  calculated  number  of  trees  available  in  the 
upstream  reach,  if  that  option  was  selected. 

3)  The  probability  results.  If  the  average  tree  height  (Ht)  is  greater  than  the  distance  between 
the  bridge  piers  (Is)  the  probability  of  at  least  one  tree  becoming  trapped  is  assumed  to  be  I  If 
Ht  <  Is  the  probability  of  each  tree  in  the  upstream  reach  successively  becoming  trapped  is 
displayed  up  to  the  total,  Nt  trees 

4)  Bridge  pier  scour  calculations.  These  values  are  calculated  assuming  that  all,  NT,  trees  are 
caught  to  give  a  worst  case  scenario  The  following  values  arc  then  displayed: 

i)  The  probability  of  NT  trees  becoming  caught 

ii)  The  alculated  debris  raft  depth 

iii)  The  percentage  of  the  channel  cross-section  blocked  by  debris. 

iv)  Bed  scour  due  to  the  pier  alone  (using  2  4  x  D) 

v)  Bed  scour  due  to  the  pier  and  debris  accumulation 

vi)  The  hydrostatic  pressure  force  on  the  pier  due  to  the  debris  accumulation 


ion 
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Test  Cases 

These  examples  are  provided  to  ensure  '>at  the  user  is  familiar  with  the  procedures  necessary 
to  run  the  programs,  and  1  '  check  that  tne  supplied  program  is  working  correctly, 

1 )  Start  the  program  by  typing  “dprob”  from  the  LWD  directory 

2)  Enter  an  output  file  name,  e  g.  “results  txt”. 

3)  Enter  the  following  values 


run  l 

run  2 

run  3 

run  4 

tree  trunk  diameter 

0.5 

0.5 

0.5 

0.5 

max  tree  diameter 

2 

2 

2 

2 

average  tree  height 

6 

6 

6 

6 

tree  estimation  option 

2 

2 

1 

1 

tree  density 

- 

0.025 

0.0125 

reach  length 

- 

- 

100 

100 

failure  width 

- 

- 

2 

4 

number  of  trees 

10 

10 

- 

- 

pier  diameter 

1 

1 

1 

1 

pier  spacing 

6 

20 

30 

40 

flow  depth 

6 

6 

6 

6 

4)  Exit  the  program  and  view  the  results  in  an  text  editor  or  word  processor. 

5)  The  results  in  runs  I  to  4  should  read  as  follows: 

Program  run  number  1 

Dl-URIS  AT  HRIlXil-  IMliR  <’ A IX/H  RATIONS 

The  tree  trunk  diameter  is  0.5  metres 
The  tree  root  wad/canopy  diameter  is  2  metres 
Thu  average  tree  height  is  ft  metres 
The  observed  number  of  trees  in  the  upstream  reach  is  10 
I  he  bridge  pier  diameter  is  I  metres 
The  length  between  bridge  piers  is  ft  metres 
Hie  flow  depth  is  6  metres 


\ 


Bridge  pier  scour  output  : 

lire  following  values  have  been  calculated  liv  assuming  all  10  trees  are  caught  at  the  bridge 
I  he  probability  of  at  least  I  out  of  10  trees  Incoming  trapped  is  100% 
lire  debris  i rill  depth  is  0.00  metres 

flic  percentage  of  the  channel  cross-section  blocked  is  10.51  % 
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Scour  calculation  arc  don  veil  from  a  model  developed  bv  Melville  and  Dongol  ( 1  ,)l)2) 
Hed  scour  due  to  the  pier  --  2.4  metres 

lied  scour  due  to  die  pier  and  debris  accumulation  -  3.4.1  metres 
Tile  pressure  force  normal  to  the  bridge  pier  is  4  Slr  Hl.i  N/in  unit  width 


Program  run  number  2 

I fl'HKIS  AT  I1KII X  IK  Pll-R  CAI  (T II  A I  IONS 

Idle  tree  trunk  diametei  is  0  5  metres 

Idle  tree  root  wad/canopy  diameter  is  2  metres 

Idle  uveruge  tree  height  is  6  metres 

Idle  ohserved  number  of  trees  in  the  upstream  reach  is  10 

1  he  bridge  pier  diameter  is  I  incUa  - 

The  length  between  bridge  piers  is  20  metres 

Idle  How  depth  is  (>  metres 


Prolmbility  Results 

Idle  probability  of  ut  least  one  out  of  10  trees  becoming  trapped  is  0.33 

Idle  probability  of  2  trees  being  trapped  is  0  1033 

Idle  probability  of  3  trees  t>eing  trapped  is  0.0304 

Idle  probability  of  4  trees  being  trapped  is  0.0083 

Idle  probability  of  5  trees  being  trapped  is  0.0021 

(lie  probability  oft-  trees  being  trapped  is  ().(XX)5 

Idle  probability  of  7  trees  being  trupped  is  H.(>762e-05 

flic  probability  of  8  trees  being  trapped  is  1  .33380-05 

Idle  probability  of  7  trees  being  trapped  is  1  482(ic-(Ki 

[die  probability  of  10  trees  being  trapped  is  8.7977c-OK 


Bridge  pier  scour  output : 

Idle  following  values  have  been  calculated  by  assuming  all  1(1  trees  are  caught  at  the  bridge 
Idle  probability  of  10  trees  becoming  trapped  is  8 7777e-0X 
Hie  debris  raid  depth  is  3  %  metres 

The  percentage  of  the  channel  cross-section  blocked  is  1 3.08  % 

Scour  calculation  are  derived  from  a  model  developed  by  Melville  and  Dongol  (1772) 

Hed  scour  due  to  the  pier  -  2.4  metres 

Hed  scour  due  to  the  pier  and  debris  accumulation  o  52  metres 
Hie  pressure  force  normal  to  the  bridge  pier  is  7.7e+04  N/ni  unit  width 


Program  run  number  J 

DIIHRIS  AT  HRIIXib.  PIHK  I'AIX’DLAd  lONS 

I  he  tree  trunk  diameter  is  0.5  metres  £ 

Idle  tree  root  wad/auiopy  diameter  is  2  metres 
Idle  average  tree  height  is  6  metres 

Idle  average  tree  density  in  the  upstream  reach  is  0.025  trees  |s-i  square  metre 
Idle  upstream  leach  length  is  100  metres 
Idle  channel  failure  width  is  2  m.  tres 

llie  calculated  maximum  mini  ei  ol  trees  in  the  upstream  reach  is  10 

Idle  bridge  pier  diameter  is  1  metres  • 

Hi  ■  length  between  bridge  piers  is  30  metres 
Idle  llow  depth  is  6  inches 

Probability  Results 

Idle  probability  of  ut  least  one  out  .4'  1 0  hccs  lecoiiuug  liapped  is  0.2  t27 
Idle  piobability  of  2  trees  leiug  liap|>ed  is  0.05tXi 
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lbe  probability  of  .3  trees  lieing  trapped  is  0.0102 
lbe  probability  of  4  trees  being  trappoi  is  O.OOl1) 
lbe  probability  of  5  trees  being  trap|ied  is  0.0003 
The  probability  of  6  trees  being  trapped  is  4.2X7')e-05 
lbe  probability  of  7  trees  being  trapped  is  5.08 1  le-0(> 
lbe  probability  of  8  trees  being  trapped  is  4.728  3e-07 
lbe  probability  of1)  trees  being  trapped  is  3.077Xe-0X 
lbe  probability  of  10  trees  being  trapped  is  I  0534e4)<> 


Bridge  pier  scour  output  : 

lbe  following  values  huve  been  calculated  hv  assuming  all  10  trees  are  caught  a!  the  bridge 
lbe  probability  of  10  trees  becoming  trapped  is  1  .0434cm)') 

Hie  debris  rut!  depth  is  3  %  meties 

lbe  percentage  of  the  cluuinel  cross-section  blocked  is  8.72  % 

Scour  calculation  are  derived  from  a  model  developed  by  Melville  and  Dongol  ( l‘)‘)2) 

Bed  scour  due  to  the  pier  =  2.4  metres 

Bed  scour  due  to  the  pier  and  debris  accumulation  ~  (>.52  metres 
lbe  pressure  force  normal  to  die  bridge  pier  is  7  7eH)4  N/m  unit  width 


Program  run  number  4 

DliBRIS  AT  HRllXii;  I'll  It  CA1.CIII  Al  lONS 

Hie  tree  trunk  diameter  is  (I  5  metres 

Hie  tree  root  wud/canopy  diameter  is  2  metres 

lbe  average  tree  height  is  (>  metres 

Hie  average  tree  density  in  the  upstream  reach  is  0.0125  bees  per  stpiare  metre 
lbe  upstream  reach  length  is  l(M)  metres 
lbe  channel  failure  width  is  4  metres 

lbe  calculated  maximum  number  of  trees  in  die  upslieam  reach  is  Id 

lbe  bridge  pier  diameter  is  I  metres 

lbe  length  between  bridge  piers  is  40  metres 

Hie  How  depdi  is  6  metres 


Probability  Results 

lbe  probability  of  ut  lens1  one  out  of  10  trees  becoming  trapped  is  0. 17‘)X 

lbe  probability  of  2  trees  being  trapped  is  0.021)'* 

lbe  probability  of  3  trees  being  Irapped  is  0.004(> 

l  be  probability  of  4  trees  lieing  trapped  is  0.0000 

lbe  probability  of  5  trees  being  trapped  is  7.5fi(i.3c-05 

flic-  probability  of  (>  trees  lieing  trapped  is  7.88  1 7e-0(i 

lbe  probability  of  7  trees  being  trapped  is  0.782  3c -07 

lbe  probability  of  8  trees  lieing  trapped  is  4.524 5e-0K 

lbe  probability  of1)  trees  lieing  It  appal  is  2.082le-0l) 

lbe  probability  of  10  trees  being  trapped  is  4.%20e-l  I 


Bridge  pier  scour  output  : 

The  following  values  have  lieen  calculated  by  assuming  all  10  trees  are  caught  at  the  luidgc 
lbe  probability  of  10  tiees  t recoining  t nipped  is  4  ‘Hi2(ie-I  I 
lbe  debris  rail  depth  is  3.90  metres 

l  be  percentage  of  the  channel  cross- section  blocked  is  0.54  % 

Scour  calculation  me  denved  from  a  iiuhIcI  develo|ied  by  Melville  and  Dongol  ( l')')2) 

Bed  seoui  doe  to  the  pin  7  4  metres 

Bed  scour  due  to  the  pici  and  dels  is  aceuiiiiilalion  0.42  inches 
lbe  op  me  force  uoi mat  to  the  In  idge  piei  is  7  7c  1 1)4  N/m  unit  width 
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Results  Discussion 

The  probability  of  at  least  one  tree  snagging  at  the  pier  in  run  I  is  assumed  to  be  100% 
because  the  tree  height  (Ht)  is  equal  to  the  pier  spacing  (Is),  while  probabilities  for  each  of  the 
ten  trees  arriving  at  the  span  in  runs  2,  and  4  are  calculated  because  Ht  is  less  than  Is 
The  probabilities  for  the  10  trees  in  runs  2,  b  and  4  can  be  plotted  as  shown  in  figure  A1  This 
plot  shows  how  the  probability  of  capture  is  dependant  upon  the  maximum  tree  diameter  to 
pier  spacing  ratio  and  demonstrates  that  the  fail  in  probability  of  capture  becomes 
asymptotically  smaller  as  the  clear  water  area  between  piers  is  reduced.  In  this  instance  the 
probabilities  all  tend  towards  zero  because  the  difference  between  pier  spacing  and  tree 
diameter  is  large. 

Figure  A1  Tree  trapping  probability  plot 


-a-  tun  2  •  tun  3  run  4 
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